AD-A16G  7*7 
UNCLASSIFIED 


ESTIMATING  RADIATED  THERMAL  ENVIRONMENTS  FROM  iAj 

ALUHINUM-LOX  (LIQUID  OXVGEN.  .  (U)  SCIENCE  AND 
ENGINEERING  ASSOCIATES  INC  SALEH  HA  B  S  CHAMBERS 
IS  DEC  84  SEA-62-GG-Q : Bi  DNA-TR-83-2  F/G  28/12  NL 


AD-A168  707 


DNA-TR-85-3 


ESTIMATING  RADIATED  THERMAL  ENVIRONMENTS  FROM 
ALUMINUM-LOX  THERMAL  RADIATION  SIMULATORS 


Burton  S.  Chambers  III 
Science  and  Engineering  Associates,  Inc. 
76  Lafayette  Street 
Salem,  MA  01970 


15  December  1984 


Technical  Report 


CONTRACT  No.  DNA  001-81-C-0020 


Approved  for  public  release; 
distribution  is  unlimited. 


TH'S  WORK  WAS  sponsored  by  the  defense  nuclear  agency 
UNDER  Rf'TAE  RMSS  CODE  B34S08  1466  G54AAXYX00001  H2b90D 


Prepared  for 
Director 

DEFENSE  NUCLEAR  AGENCY 
Washington,  DC  20305- 1000 


UNCLASSIFIED _ 

security  classification  of  this  page 


REPORT  DOCUMENTATION  PAGE 


la  REPORT  SECURITY  CLASSIEICATIOlM 

UNCLASSIFIED 


2a  SECURITY  CLASSIFICATION  authority 


2b  DECLASSIFICATION /DOWNGRADING  SCHEDULE 

N/A  since  UNCLASSIFIED 


a  performing  organization  report  number(S) 


SEA  62 -00-Q : 0 1 


6c  ADDRESS  (G(y,  State,  and  ZIP  Code) 

76  Lafayette  Street 
Salem,  MA  01970 


6b  OFFICE  SYMBOL 
(If  applicable) 


Form  Approved 
OMBNo.  0704 -0188 
Exp.  Date:  lunZO.  1986 


1b  RESTRICTIVE  MARKINGS 


3  distribution /availability  of  report 

Approved  for  public  release;  distribution 
is  unlimited. 


5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

DNA-TR-85-3 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Director 

Defense  Nuclear  Agency 


7b  ADDRESS  (C/fy,  State,  and  ZIP  Code) 

Washington,  DC  20305-1000 


8a  NAME  OF  FUNDING  -  SPONSORING 
ORGANIZATION 


8c  ADORESS  (C/ty,  State,  and  ZIP  Code) 


8b.  OFFICE  SYMBOL  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable) 

DNA  001-81-C-0020 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM 
ELEMENT  NO 

62715H 


PROJECT 

TASK 

NO 

NO 

G54AAXY 

X 

WORK  UNIT 
ACCESSION  NO 

DH005106 


11  TITLE  (Include  Security  Clauification) 

ESTIMATING  RADIATED  THERMAL  ENVIRONMENTS  FROM  ALUMINUM-L0X  THERMAL  RADIATION  SIMULATORS 


12  PERSONAL  AUTHOR(S) 

Burton  S.  Chambers  III 


13a  type  OF  REPORT 

Technical 


16  supplementary  notation 

This  work  was  sponsored  by  the  Defense  Nuclear  Agency  under  RDT&E  RMSS  Code 
B345081466  G54AAXYX00001  H2590D. 


14  DATE  OF  REPORT  (Year.  Month.  Day)  15.  PAGE  COUNT 

1984,  December  15  190 


FIELD 

GROUP 

14 

2 

20 

13 

COSATI  CODES 


SUB-GROUP 


18  SUBJECT  TERMS  {Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Nuclear  Weapons  Effects  TRS  Free-Field  Data 

Thermal  Radiation  Simulators(TRS)  MILL  RACE 

Simulation  of  Thermal  Radiation  Airblast/TRS  Interference 


19  ABSTRACT  {Continue  on  reverse  if  necessary  and  identify  by  block  number ) 

This  report  describes  the  early  modeling  work  done  to  estimate  the  radiated  thermal  environ¬ 
ment  generated  by  Aluminum-LOX  thermal  radiation  simulators.  The  ability  to  calculate  the 
environment  produced  by  these  simulators  is  necessary  to  effectively  plan  combined  airblast/ 
thermal  simulations  of  nuclear  weapon  effects.  Included  are:  (1)  an  overview  of  the  model¬ 
ing  effort  and  how  a  model  is  calibrated,  (2)  comparisons  with  MILL  RACE  results,  (3)  the 
models,  and  (4)  preliminary  work  on  airblast/TRS  interference.  The  primary  result  of  the 
effort  was  the  development  of  a  model  that  potentially  can  accurately  predict  the  environ¬ 
ment  from  a  TRS.  It  was  also  concluded  that  much  more  effort  must  go  into  making  repeatable 
and  reliable  TRS  simulators,  and  a  significant  calibration  program  is  then  warranted. 


A8STRAC'  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


20  DISTOIBUT. ON  AVAILABILITY  OF  ABSTRACT 

L3  UNCLASSiFiED'UNL.MlTED  □  SAME  AS  RPT  □  one  USERS 


22a  NAME  OF  RESPONSIBLE  NDiVIDUAl  22b  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOL 

Betty  L.  Fox  (202)  325-7042  DNA/STTI 


DO  FORM  1473,  94  MAR  83  APR  edition  may  be  used  until  exhausted  SECURITY  CLASS1F1CAT  ON  OF  -h,s  page 

All  other  editions  are  obsolete  UNCLASSIFIED 


6a.  NAME  OF  PERFORMING  ORGANIZATION  (Continued) 

Under  purchase  order  (No.  15-840006  65)  from: 

Science  Applications  Inti  Corp. 

6c.  ADDRESS  (Continued) 

P.0.  Box  1303 
McLean,  VA  22102-1303 

18.  SUBJECT  TERMS  (Continued) 

Modeling  of  Simulated  Environments 
Airblast/Thermal  Synergistic  Effects 


SUMMARY 


The  work  reported  herein  was  conducted  in  support  of  the  use  of 
thermal  radiation  simulators  (TRS)  on  the  Defense  Nuclear  Agency  (DNA)  high 
explosive  (HE)  event  MILL  RACE.  The  type  of  TRS  used  on  MILL  RACE  was  the  so- 
called  flame  or  torch  TRS,  which  consists  of  a  set  of  nozzles  each  burning 
aluminum  powder  in  sprayed  liquid  oxygen  (LOX). 


Guidance  and  technical  assistance  were  provided  to  DNA  Field  Command 
and  in  turn  to  any  MILL  RACE  experimenters,  who  would  be  utilizing  the  TRS. 
Models  of  the  thermal  radiation  were  developed  and  then  estimates  were  made  of 
the  expected  thermal  environments  at  various  locations.  New  TRS  data  were 
incorporated  into  a  data  base  which  was  then  used  to  improve  the  mode1 
predictions.  Recommendations  were  made  indicating  where  calorimeters  should 
be  placed  in  order  to  best  support  the  modeling  effort.  Unfortunately,  other 
requirements,  such  as  the  determination  of  proper  operating  conditions  for  the 
TRS  nozzles,  generally  precluded  use  of  the  recommendations  concerning  data 
collection  for  improving  the  models. 


In  addition  to  the  estimation  of  thermal  environments,  this  effort 
contained  a  task  to  investigate  interference  to  the  HE-generated  air  blast 
from  thermal  effects  generated  by  the  TRS  units.  Of  primary  concern  were 
effects  that  would  adversely  affect  planned  air  blast  loading.  If  the 
interference  effects  were  considered  to  be  significant,  suggested  alternative 
approaches  were  to  be  made  toward  meeting  the  objectives  of  the  combined 
thermal-air  blast  experiments. 


Finally,  the  models  were  again  refined  as  a  part  of  this  effort, 
after  the  MILL  RACE  event  occurred.  The  objective  of  this  latter  task  was  to 
improve  the  model  developed  earlier  by  incorporating  two  modifications:  (I) 
represent  each  flame  as  an  inverted  cone  rather  than  a  right  circular 
cylinder,  and  (2)  account  for  the  change  in  view  angle  resulting  from  flame 
obscuration. 
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This  effort  produced  a  number  of  results  and  conclusions  during  the 
course  of  estimating  the  radiated  thermal  environments  from  aluminum-liquid 
oxygen  (LOX)  thermal  radiation  simulators  (TRS).  Most  importantly,  it 
produced  models  of  TRS  generated  peak  flux  that  represent  those  same 
environments  in  the  field.  The  inverted  cone  model  is  the  best  developed  to 
date  when  flame  obscuration  and  an  accurate  treatment  of  view  angle  are 
included. 


It  is  now  possible  to  use  one  calibration  of  the  inverted  cone  model 
to  represent  all  nozzle  arrangements  used  thus  far.  Consequently,  the  data 
base  may  now  be  of  sufficient  size  to  allow  certain  statistical  inferences  to 
be  made  about  the  inherent  variability  of  using  the  TRS  in  the  field. 

In  the  area  of  TRS-Air  Blast  interference,  the  results  are  less 
clear.  The  effort  conducted  is  largely  incomplete,  and  additional  work  is 
required  to  resolve  the  issues.  Calculations  indicate  that  incipient  blowoff 
of  ground  material  was  probable  on  MILL  RACE.  Other  calculations  have  also 
indicated  that  incipient  ground  blowoff  is  probable  at  the  DIMA  TRS  site  at 
Kirtland  Air  Force  Base  (KAFB)  in  Albuquerque,  New  Mexico.  What  looks  like 
the  "steam  porch"  phase  of  ground  blowoff  has  been  seen  in  the  photographs 
taken  at  the  TRS  site.  The  effect  of  any  resultant  hot  thermal  layer 
immediately  above  the  ground  on  air  blast  propagation  was  not  investigated. 

It  was  also  concluded  that  modeling  a  TRS  flame  as  merely  a  hot  clean 
flame,  i.e.  one  containing  no  particulates,  will  not  accurately  simulate 
results  from  MISERS  BLUFF.  During  MISERS  BLUFF  the  overpressure  was  reduced 
on  one  of  the  targets  placed  immediately  behind  the  TRS  units.  Numerical 
calculations  performed  for  overpressures  up  to  10  psi  and  for  flame  widths  up 
to  two  meters  wide  demonstrated  that  other  effects  will  need  to  be  considered. 
Additional  work  is  recommended  to  quantitatively  understand  the  influence  of 
particulates  on  the  effect  seen  during  MISERS  BLUFF. 

Our  recommendations  in  the  area  of  TRS  modeling  are  that  much 
experimental  work  and  some  modeling  work  need  to  be  done.  In  the  area  of  TRS- 
airblast  interference,  our  recommendation  is  to  determine  how  well  the  current 
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interference  mitigation  techniques  will  find  acceptance  in  the  nuclear 
survivability  and  vulnerability  community.  The  current  technique  is  to  fire 
the  TRS  early  enough  to  allow  the  flame  and  its  combustion  products  to  rise 
above  the  target  so  that  they  will  not  interfere  with  the  air  blast.  For  some 
nuclear  scenarios,  this  approach  will  introduce  a  significant  cooling  period 
before  air  blast  arrival. 

A  substantially  more  exhaustive  and  careful  effort  to  obtain 
calorimeter  data  needs  to  be  performed.  This  recommendation  is  made  with  the 
knowledge  that  the  data  base  essentially  doubled  during  September  1982  when 
DNA  and  SAI  conducted  a  Calibration  Program  of  the  DNA  TRS  site  at  Kirtland 
AFB.  Even  with  those  data,  significant  additional  data  obtained  under  more 
controlled  conditions  is  necessary  if  the  community  is  to  be  expected  to  place 
great  reliance  on  modeled  results. 

More  thorough  data  interpretation  and  analysis  is  recommended.  At 
least  four  individuals  have  read  the  "peak"  flux  off  oscilliscope  traces 
that  have  been  entered  into  the  DNA  data  base.  Now  that  these  data  are  being 
recorded  digitally,  the  detailed  temporal  waveforms  can  be  made  available  for 
meaningful  data  analysis.  There  has  been  no  consistent  approach  applied  to 
evaluating  these  data.  In  some  cases  data  were  withheld  as  not  being 
representative. 

The  TRS  development  to  date  has  been  geared  toward  improving  the 
flame  generation.  Little  emphasis  has  been  given  to  evaluating  reliability  or 
repeatability.  We  strongly  recommend  that  one  design  be  fixed  for  future  use 
in  the  DNA  HE  program,  and  that  that  design  be  calibrated  under  realistic 
field  conditions. 

Finally,  we  regretably  observe  that  much  difficulty  in  accessing  the 
DNA-provided  computational  facilities  at  KAFB  made  it  virtually  impossible  to 
continue  the  investigation  of  TRS-air  blast  interference.  It  is  recommended 
that  DNA,  not  a  contractor,  more  closely  monitor  difficulties  that  contractors 
are  having  when  using  government  furnished  equipment.  Whether  this  would  lead 
to  any  improvement  in  throughput  is  debatable;  nevertheless,  we  feel  it  is 


important  to  improve  this  communication  of  the  problems  we  have  been  having, 
so  that  DNA  management  might  have  representative  information  available  when 
decisions  are  made. 
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SECTION  I 


INTRODUCTION 

The  effort  reported  herein  was  conducted  in  support  of  TRS 
experiments  on  MILL  RACE  and  consisted  of  a  pretest  and  posttest  phase.  The 
objectives  for  the  pretest  phase  were:  (I)  to  provide  guidance  and  technical 
assistance  to  DNA  Field  Command  and  the  MILL  RACE  TRS  experimenters,  and  (2) 
investigate  possible  interference  of  the  TRS  units  with  the  MILL  RACE  air 
blast,  and  then  to  suggest  alternative  approaches  for  meeting  experimental 
objectives,  if  interference  effects  were  expected  to  be  significant.  The 
objective  of  the  posttest  phase  was  to  improve  the  modeling  of  the  TRS 
radiated  environment,  where  flames  partially  obscure  one  another,  by  modifying 
the  earlier  obscuration  model  as  follows:  (I)  represent  the  flame  as  an 

inverted  cone  rather  than  a  right  cylinder,  and  (2)  account  for  the  change  in 
view  angle  resulting  from  flame  obscuration. 


I- 1  THE  DNA  THERMAL  SIMULATION  PROGRAM. 

An  important  difference  between  explosions  of  nuclear  weapons  and 

chemical  high  explosives  (HE)  is  the  fraction  of  energy  released  in  the 
electromagnetic  spectrum.  Since  a  nuclear  explosion  starts  at  a  much  higher 
energy  state,  the  transport,  absorption,  and  reradiation  of  this  energy 
becomes  quite  significant  (Reference  I).  The  chemical  HE,  on  the  other  hand, 
releases  most  of  its  energy  in  mechanical  form,  ie.  air  blast. 

During  the  nuclear  explosion  a  large  amount  of  energy  is  reradiated 

from  the  X-ray  generated  fireball  in  the  form  of  light  and  heat;  both  forms 
are  commonly  referred  to  as  thermal  radiation.  This  thermal  radiation, 
generally  arriving  before  the  shock  wave,  can  precondition  an  exposed  weapon 

system,  so  that  its  response  to  air  blast  may  be  modified  from  what  it  would 
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be  to  air  blast  alone. 


The  combination  of  air  blast  and  thermal  radiotion  loading  can  cause 
synergistic  effects  (ones  more  severe  than  what  would  be  expected  from  both 
loadings  operating  independently).  These  synergistic  effects  have  been 
observed  on  various  military  equipment  (Reference  2).  Until  fairly  recently, 
however,  these  effects  have  been  impossible  to  generate  in  the  field,  without 
above-ground  nuclear  testing,  due  to  the  lack  of  thermal  simulators  that 
produce  significant  radiative  output  over  large  areas. 

DNA  has  been  providing  simulated  nuclear  weapon  environments  with  HE 
events  to  investigate  various  nuclear  weapon-related  effects.  The  simulated 
environments  provided  hove  been  air  blast  and  ground  motion.  The  HE  events 
have  been  conducted  for  many  years  to  support  the  United  States  military 
services,  and  in  some  cases,  our  allies.  More  recently,  the  Shock  Physics  Test 
Division  (SPTD)  at  DNA  has  been  supporting  research  that  is  identifying 
various  ways  to  simulate  intense  thermal  radiation,  characteristic  of  nuclear 
weapons,  in  order  to  study  both  synergistic  thermal-air  blast  effects  and 
thermal-only  effects. 

DNA  is  sponsoring  both  the  development  of  large  area  thermal 
irradiation  simulators  to  allow  observation  of  full  scale  effects  and  the 
development  of  smaller  laboratory  thermal  simulators  that  provide  better 
simulation  fidelity.  Effective  use  of  both  types  of  these  facilities  should 
allow  significant  advances  in  our  knowledge  of  nuclear  weapon  thermal  effects 
and  synergistic  thermal-air  blast  effects. 

The  first  significant  large  area  irradiation  system,  initiated  seven 
years  ago,  was  the  so-called  bag  system  developed  by  Dishon  and  Lottery  of 
Science  Applications,  Inc.  (SAD  for  DNA.  The  bag  system  consisted  of  an 
ensemble  of  mylar  bags,  each  bag  being  a  cylinder,  bemisphericaily  capped  at 
both  ends,  about  6  meters  high  with  a  diameter  of  1.5  meters.  Each  bag  was 
filled  with  pure  gaseous  oxygen  into  which  was  sprayed  aluminum  powder.  This 
mixture  was  then  ignited  with  a  device  resembling  a  common  railroad  flare. 


The  bag  system  was  fielded  during  the  late  seventies  at  a  site  used 

by  DNA  on  Kirtland  AFB.  Various  targets  of  military  sig  sificance  were  exposed 
to  thermal-only  environments.  These  tests  allowed  qualitative  assessments  of 

the  effects  induced  from  large  area  irradiation  onto  targets  of  significant 

size.  Effects  such  as  aircraft  panel  buckling,  which  would  be  hard  to  obtain 

on  small  tests,  were  observed.  This  type  of  data  was  felt  to  be  useful  for 
validating  large  complex  structural  analysis  computer  codes. 

One  ensemble  of  TRS  bags  was  fielded  on  DNA  HE  event  MISERS  BLUFF 
where  it  became  obvious  that  the  use  of  the  TRS  actually  interferred  with  the 
air  blast  loads  delivered  to  the  exposed  U.S.  Army  S-280  hardened  shelter. 
The  bag  system  had  been  placed  in  close  proximity  to  the  shelter,  between  it 
and  the  HE.  When  the  air  blast  propagated  through  the  hot  combustion  products 
of  the  bags,  its  peak  overpressure  was  significantly  reduced.  This  experiment 
provided  some  of  the  motivation  for  seeking  an  improved  way  of  burning  the 
aluminum  powder. 

A  significant  advancement  in  simulators  was  achieved  when  Dishon  and 
his  co-workers  at  SAI  successfully  developed  the  flame  TRS  system  under  DNA 
sponsorship.  The  flame  system  is  now  recognized  as  a  significant  improvement 
over  the  bag  system,  not  surprisingly  since  it  was  designed  to  avoid  many  of 
the  undesirable  features  of  its  predecessor. 

The  flame  system  consists  of  a  set  of  nozzles,  vaguely  similar  to 
rocket  engines,  but  directed  upwards,  that  burn  aluminum  powder  in  an  oxygen 
rich  environment.  The  aluminum  powder  is  fluidized  in  nitrogen  gas  to 
facilitate  flow  through  the  plumbing  from  supply  to  nozzle.  Liquid  oxygen 
(LOX)  is  also  sprayed  into  this  mixture  in  the  nozzle.  The  resulting 
combustible  mixture  is  ignited  with  a  flame  that  is  burning  outside  the 
nozzle.  This  outside  flame  is  lighted  before  the  aluminum  and  liquid  oxygen 
mixture  is  flowing,  and  its  fuel  is  usually  propane.  After  the  aluminum  and 
oxygen  mixture  is  ignited,  a  bright  flame  is  produced  that  can  burn  until  the 
materials  being  injected  are  exhausted. 


Due  to  the  time  taken  for  the  burning  of  the  aluminum  particles,  the 


initial  momentum  associated  with  the  spray,  and  buoyancy  forces,  the  flame 
from  each  nozzle  typically  reaches  a  significant  height  of  three  to  eight 
meters  depending  on  how  rapidly  aluminum  is  introduced  into  the  system.  The 
radiating  part  of  each  flame,  for  a  5  kilogram  per  second  burn  of  aluminum, 
appears  approximately  as  an  inverted  cone  of  height  about  6  meters  and 
diameter  of  5  meters  at  the  cone's  "base"  (the  nozzle  is  at  the  cone's  apex). 

This  latter  simulator,  the  Ai-LOX  flame  system,  was  fielded  on  MILL 
RACE,  and  is  the  system  that  has  been  modeled  during  this  effort. 

1-2  THE  DNA  THERMAL  SIMULATION  MODELING  PROGRAM. 

DNA  has  also  been  sponsoring  research  on  calibrating  these  AI-LOX 
flame  sources  to  allow  their  use  in  the  field,  either  as  stand-alone  systems 
or  in  conjunction  with  large  HE  events.  The  experimental  calibration  programs 
have  been  complemented  with  modeling  efforts  to  allow  interpolation,  and  in 
some  cases  extrapolation,  of  the  environments. 

In  preparation  for  the  MILL  RACE  experiment,  the  effort  reported 
herein  was  initiated  to  provide  guidance  and  technical  assistance  to  DNA  Field 
Command  and  to  the  MILL  RACE  experimenters  who  would  be  utilizing  the  AI-LOX 
flame  source.  Specifically,  estimates  were  made  of  the  thermal  environments 
expected  from  the  TRS  units  at  various  experiment  locations  and  then  provided 
to  FCDNA  for  planning  the  use  of  the  TRS  units. 

The  estimates  were  based  on  early  models  of  the  bag  system,  that  had 
been  developed  from  a  very  limited  data  base.  As  data  generated  under  another 
DNA  contract  were  made  available,  they  were  incorporated  into  the  DNA  TRS  data 
base  and  were  used  to  expand  and  improve  the  existing  models. 

When  the  test  bed  configuration  was  altered,  or  when  significantly 
different  (and  better)  estimates  of  the  AI-LOX  TRS  environment  became  available, 
updated  estimates  were  provided  to  FCDNA  and  these  were  in  turn  given 


to  the  MILL  RACE  experimenters.  Recommendations  were  also  made,  on 
a  limited  basis,  where  calorimeters  should  be  placed  in  order  to  improve  the 
modeling  effort.  Other  requirements,  however,  generally  precluded  use  of 
these  recommendations  by  the  TRS  developers. 

1-3  THE  DNA  THERM AL-AIR  BLAST  INTERFERENCE  PROGRAM. 

In  addition  to  the  estimation  of  thermal  environments,  this  effort 
included  a  task  to  investigate  possible  modifications  to  (or  interference 
with)  the  HE-generated  air  blast  by  thermally  induced  effects  generated  by  the 
TRS  units.  Of  primary  concern  were  adverse  effects  on  planned  air  blast 
loading,  as  was  seen  on  MISERS  BLUFF  (Reference  2).  If  the  interference 
effects  were  considered  significant,  suggested  alternative  approaches  for 
meeting  the  objectives  of  the  combined  thermal-air  blast  experiments  were  to 
be  made. 

It  is  necessary  to  place  the  flames  relatively  close  to  the  targets 
to  be  irradiated  since  the  thermal  output  of  even  these  powerful  TRS 
simulators  is  small  in  comparison  to  that  of  a  nuclear  burst  (approximately 
five  orders  of  magnitude  below  a  kiloton  of  energy).  This  arrangement 
produces  combined  air  blast  and  thermal  effects  under  certain  conditions. 
These  effects  can  approximate  those  that  would  be  seen  during  an  actual 
nuclear  explosion.  However,  if  the  TRS  simulator  is  not  sufficiently  strong 
to  allow  its  placement  far  enough  away  from  the  irradiated  target,  it  is 
possible  that  the  combined  effects  generated  may  only  fortuitously  simulate 
real  ones. 

An  example  was  the  experience  on  MISERS  BLUFF  where  an  S-280  U.S.  Army 
Shelter  was  exposed  to  both  thermal  and  air  blast  (Reference  2).  The  only 
available  large  scale  simulator  for  MISERS  BLUFF  was  the  bag  TRS.  It  had  to 
be  placed  very  close  to  the  target  in  order  to  be  able  to  deliver  the  required 
environment.  Because  the  timing  between  the  onset  of  thermal  irradiation  and 
air  blast  arrival  was  set  to  simulate  a  particular  nuclear  scenario  with  this 
MISERS  BLUFF  experiment,  the  air  blast  arrived  while  the  TRS-generated 


fireball  was  still  in  front  of  the  shelter.  When  the  air  blast  passed  through 
this  hot  region,  which  was  essentially  touching  the  target  surface,  the 
received  air  blast  overpressure  was  appreciably  reduced.  This  effect  was 
considered  undesirable  by  the  U.S.  Army. 

On  MILL  RACE,  however,  an  attempt  was  made  to  mitigate  the  direct 
influence  of  the  flame  on  the  air  blast  peak  overpressure  by  allowing 
sufficient  time  for  the  flame  to  just  rise  out  of  the  way  before  air  blast 

arrival.  The  length  of  time  necessary  for  this  to  occur  was  determined  during 

a  series  of  small  yield  HE  shots  at  the  TRS  test  site  at  Kirtland  AFB.  This 
modest  HE  test  program,  called  the  Pentolite  series,  also  determined  if  the 
air  blast  would  be  modified  by  passage  through  the  region  of  space  previously 
occupied  by  the  flame.  This  approach  of  avoiding  thermal  "interference"  with 
air  blast  takes  advantage  of  the  expected  time  delay  in  the  real  world  case 
and  is  the  current  planned  approach  for  use  of  TRS  units  during  future  HE 

tests. 

However,  careful  consideration  of  the  consequences  of  this  approach 
are  advisable  since  its  use  could  lead  to  other  problems  of  interpretation. 

For  example,  if  the  target  cools  too  much  before  air  blast  arrival,  the  test 
may  not  meaningfully  represent  the  worst  possible  effects.  One  might 
incorrectly  conclude  that  some  system  is  less  vulnerable  to  combined  effects 

than  it  really  is. 

1-4  TRS  UNITS  FIELDED  AT  MILL  RACE. 

Four  flame  systems  were  eventually  fielded  for  the  DNA  HE  event,  MILL 
RACE  (Reference  3).  Throughout  this  report,  the  systems  are  referred  to  as: 
(I)  UK- 1,  (2)  UK-2,  (3)  BRL,  and  (4)  Navy.  Each  of  the  four  TRS  systems 

consisted  of  one  line  of  equally  spaced  nozzles.  The  UK- 1  system,  fielded  for 
the  use  of  the  United  Kingdom  experimenters,  consisted  of  eight  nozzles  of 
seven  feet  separation.  The  UK-2  system,  also  for  the  United  Kingdom, 

consisted  of  four  nozzles  spaced  at  the  same  seven  foot  interval.  The  BRL 
system,  fielded  primarily  for  use  by  the  U.  S.  Army  Hardened  Army  Tactical 


Shelter  (HATS)  program,  consisted  of  four  nozzles  spaced  three  and  one-half 
feet  apart.  Finally,  the  Navy  system,  fielded  for  various  U.S.  Navy 
experiments,  consisted  of  four  nozzles  spaced  ten  and  one-half  feet  apart. 

1-5  REPORT  OUTLINE. 

The  principal  product  of  this  effort  has  been  the  development  of 
models  that  describe  the  radiated  thermal  environment  in  the  field  from  AI-LOX 
TRS  systems  burning  at  the  rate  of  5  kilograms  of  aluminum  per  second. 
Although  numerous  calculations  were  performed  that  were  useful  in  helping  DNA 
site  the  TRS  experiments  on  MILL  RACE,  none  of  these  are  difficult  nor 
expensive  to  replicate,  and  therefore,  there  is  little  reason  to  document  them 
beyond  the  numerous  technical  memoranda  that  have  already  been  generated.  It 
is  worthwhile,  however,  to  present  what  has  been  learned  during  the  modeling 
effort,  since  this  will  have  lasting  value  to  those  who  will  use  the  TRS  on 
future  HE  tests  and  in  large  shock  tubes. 

Therefore,  the  bulk  of  the  report  discusses  the  models  of  TRS 
generated  environments.  Section  2  presents  an  overview  of  the  modeling  effort 
and  discusses  how  a  model  is  calibrated  before  being  used.  Section  3  presents 
comparisons  of  results  with  MILL  RACE.  Section  4  presents  the  models  and 
algorithms  that  were  developed.  Section  5  summarizes  the  work  done  on  TRS 
interference  with  air  blast.  The  results  given  are  preliminary  in  nature, 
since  not  all  of  the  issues  have  been  resolved.  Finally,  Sections  6  and  7 
present  our  conclusions  and  recommendations. 


SECTION  2 


MODELING  TRS  RADIATED  ENVIRONMENT 

This  section  presents  an  overview  of  the  modeling  activity  and  the 
goals  for  this  effort.  A  more  detailed  presentation  of  the  actual  models  used 
to  estimate  the  thermal  environments  produced  by  the  DNA  AI-LOX  TRS  is  given 
in  Section  4. 

All  of  the  various  versions  of  the  models  discussed  here  have 
historically  been  referred  to  as  ALFGEE,  which  is  an  acroymn  for  Aluminum-LOX 
Field  Generated  Environment  Estimator.  Basically,  ALFGEE  is  simply  a 
numerical  simulation  that  represents  the  current  capability  to  estimate  TRS- 
generated  thermal  environments. 

In  the  early  attempts  to  model  TRS  radiated  environments,  the  goal 
was  to  develop  ways  of  estimating  these  data  within  some  acceptable  accuracy. 
Consistency  of  model  parameters  with  the  physical  observables  was  desirable 
but  not  essential.  Later  it  became  apparent  that  more  reliance  would  be 
placed  on  the  modeling  capability  because  of  the  costs  associated  with 
measurement,  and  therefore  the  goal  shifted  to  developing  a  more  universal 
model,  one  that  agreed  better  with  the  physical  observables. 

The  most  important  motivation  for  improving  the  models  was  due  to 
their  use  in  HE  events.  Model  predictions  are  used  to  lay  out  experiment 
locations.  Inadequacies  in  estimation  can  translate  into  experimental 
requirements  being  compromised. 


2-1  MODELING  FLUX. 


When  the  TRS  units  are  ignited  the  radiated  power  rises  rapidly  to 


its  peak  in  a  few  tenths  of  a  second.  The  power  from  the  flames  then  remains 
at  the  peak  for  a  duration  close  to  the  burn  time,  typically  a  second  or  so, 
and  then  decays  rapidly.  Usually,  the  result  is  that  most  of  the  received 
energy  comes  from  the  flames  while  they  are  radiating  at  their  peak  power. 

The  resultant  environment,  that  has  been  modeled  during  this  and 
earlier  efforts,  is  the  flux  received  at  a  surface  placed  in  the  field  facing 
the  flames.  The  flux  is  modeled  in  calories  per  square  centimeter  per  second. 
Furthermore,  this  modeling  has  been  restricted  to  times  when  the  flame  has 
stabilized  in  shape  and  power  corresponding  to  when  the  flux  is  at  or  near  its 

peak  value. 

A  suggested  alternate  approach  is  to  model  fluence  measured  in 
calories  per  square  centimeter.  However,  we  have  chosen  to  not  explicitly 
model  the  fluence.  Instead,  it  is  calculated  as  the  product  of  the  burn  time 
and  the  modeled  peak  flux.  This  only  approximates  the  fluence  since  the  flame 
dimensions  depend  on  time.  Errors  associated  with  the  chosen  approach  have 
not  been  evaluated  yet  for  a  variety  of  reasons.  Furthermore,  the  need  to 
independently  model  the  fluence  has  not  been  urgent,  since  the  TRS  units  have 
normally  been  used  to  generate  relatively  long  pulses. 

While  advantages  may  exist  to  model  the  fluence,  reasons  exist  to 
prefer  modeling  the  flux,  if  a  choice  has  to  be  made.  First,  modeling  fluence 
requires  the  time  dependency  of  all  the  parameters;  otherwise  the  burn  time 

would  have  to  be  included  in  the  model  in  a  complicated  way.  Second,  any 

smoke  generated  during  an  experiment  could  corrupt  the  late  time  flux,  and 
hence  the  fluence,  in  a  way  that  greatly  complicates  the  analysis.  Sufficient 
photographic  data  do  not  exist  to  establish  which  calorimeters  have  been 
obscured  by  smoke,  hence  past  data  may  have  limited  value.  Third,  the  data 
taken  early  in  the  program  suffered  from  a  thermal  effect  on  the  calorimeter 

test  stands.  The  poles,  holding  the  calorimeters,  would  tip  back  as  they  were 
irradiated  as  a  result  of  asymmetrical  thermal  stress.  The  diameters  of  the 
poles  used  in  the  early  tests  were  too  small.  This  latter  effect  was  studied 
and  a  summary  is  presented  as  Appendix  A. 


All  of  the  effects  listed  above  tend  to  alter  the  flux  data 

especially  at  later  times.  Therefore,  we  chose  to  include  the  early  time 

(peak  at  earliest  time)  data  in  the  modeling  data  base  expecting  that  this 

would  be  sufficient  for  our  needs.  Modeling  the  rest  of  the  effects,  that  may 
be  needed  some  day,  was  not  considered  a  reasonable  approach  for  the  initial 

phase  of  the  modeling  effort.  At  that  time  few  data  existed  anyway,  so 
little  benefit  would  have  resulted. 


2-2  INTRODUCTION  TO  THE  ALFGEE  MODELS. 

An  ALFGEE  model  consists  of  relationships  and  parameters  that 
determine  the  received  thermal  radiation  at  some  detector  given  additional 
information  about  the  number  of  nozzles,  their  spacing,  and  such.  An  example 
of  a  relationship  is  the  way  the  radiated  power  per  unit  height  varies  as  a 
function  of  height.  An  example  of  a  parameter  is  the  height  of  the  radiating 

part  of  the  flame.  A  relationship  is  typically  described  by  an  equation  and 

one  or  more  parameters.  A  parameter  is  described  by  a  single  numerical  value. 

The  calculation  of  the  received  environment  in  the  field  is 
straightforward  once  the  geometrical  shape  of  a  TRS  flame  has  been 
established,  if  its  radiance  is  known.  The  irradiance  is  calculated  at 
desired  positions  in  the  field  upon  some  imaginary  surface,  the  normal  of 

which  is  oriented  in  some  specified  direction,  by  simply  integrating  over  the 
radiating  elements,  accounting  for  the  viewing  angles  to  both  the  detector  and 
each  element.  Atmospheric  attenuation  is  insignificant  for  most  of  the  energy 
for  the  distances  involved,  and  therefore,  is  not  included. 

A  number  of  approaches  have  been  used  during  the  modeling  effort  to 
represent  the  geometrical  nature  of  the  flames  when  calculating  irradiance. 
The  first  representation  was  a  linear  array  of  point  sources  each  radiating 
with  constant  power.  Initially,  it  was  easy  to  modify  parameters  to  produce 
adequate  agreement  with  the  limited  data.  This  was  especially  true  because 

many  tests  were  conducted  with  the  same  calorimeter  layout,  which  although 
providing  a  measure  of  how  well  the  TRS  development  was  proceeding,  did  not 
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yield  information  at  other  locations  needed  to  stress  the  modeling  task.  Many 
combinations  of  parameters  can  produce  good  agreement. 


The  first  refinement,  motivated  by  the  eventual  taking  of  data  at 
other  heights,  was  to  use  the  same  linear  array  of  point  sources  as  in  the 

first  representation  but  with  a  linear  failoff  in  power  with  height.  This 
second  model  served  the  community  well  until  interest  existed  in  placing 
calorimeters  off  to  the  side  of  the  TRS  units.  It  subsequently  became  clear 

that  the  model  needed  to  account  for  the  effect  of  flames  obscuring  a 
detector's  view  of  other  flames.  A  cursory  analysis  indicated  that 
considering  this  obscuration  would  substantially  improve  agreement. 

The  radiating  right  circular  cylinder  shape  was  selected  as  the  third 

representation  for  its  simplicity  in  dealing  with  obscuration.  It  proved  to 

be  quite  acceptable  until  data  became  available  for  the  early  test  firings  of 
the  MILL  RACE  nozzles.  Before  those  events,  all  the  data  had  been  taken  for  a 
nozzle  spacing  of  3.5  feet.  With  the  early  test  firing  data  available  (two 
sets  of  7  feet  spacing  and  one  set  of  10.5  feet  spacing),  it  became  clear  that 
it  would  be  necessary  to  adjust  model  parameters  for  each  nozzle  spacing. 

Further  study,  however,  suggested  that  it  might  be  possible  to  avoid 
recalibrations  for  each  spacing  by  changing  the  representation  of  the  flame's 
geometrical  shape.  Consequently,  additional  work  was  performed  to  improve  the 
agreement  between  data  and  models  subsequent  to  the  MILL  RACE  event.  This 
additional  work  resulted  in  the  fourth  representation  that  was  designed  to 
approximate  the  actual  flame  shape.  Video  and  photographic  records  show 
flames  that  look  like  inverted  cones  with  dense  clouds  of  combustion  products 
above  them.  Therefore  the  final  representation  was  patterned  as  a  radiating 
inverted  cone  capped  with  an  opaque  right  cylinder.  The  power  per  unit  length 
relationship  was  modified  to  reflect  this  newer  geometry. 

The  use  of  the  inverted  cone  representation  has  been  found  to 
adequately  replicate  the  experimental  data.  It  appears  that  a  single 
calibration  of  the  model  provides  reasonable  agreement  with  the  data  base  for 
a  nozzle  that  burns  5  kilograms  of  aluminum  per  second  for  many  nozzle 


spacings.  Since  the  inverted  cone  model  was  the  best  one  developed  during 
this  effort,  it  is  described  in  more  detail  in  the  following  subsections. 

An  additional  point  must  be  presented  before  proceeding  to  the 
discussion  of  the  inverted  cone  model.  As  a  result  of  our  review  of  the  data, 
we  conclude  that  there  is  reason  to  believe  that  each  TRS  unit,  by  itself,  is 
reproducible  if  fired  under  the  same  operating  cor  itions  and  characteristics. 
There  is  much  variability  in  the  data,  however,  and  some  of  this  we  suspect  is 
due  to  differences  among  the  TRS  units,  for  example,  nozzle  design. 

2-3  INTRODUCTION  TO  THE  INVERTED  CONE  MODEL. 

After  MILL  RACE  the  inverted  cone  model  was  developed  and  algorithms 
for  flame  obscuration  effects  were  included.  Each  flame  was  treated  as  a 
radiating  inverted  cone  capped  by  an  opaque  cylinder,  consistent  with 
available  photographic  data.  The  radiating  portion  is  considered  to  be  a 
diffuse  source.  Capping  the  flame  with  the  opaque  cylinder  approximates  the 
influence  from  smoke  generated  by  the  flames.  In  addition,  the  effect  due  to 
change  in  view  angle  resulting  from  flames  obscuring  one  another  was 
considered  and  added  to  the  model.  The  energy  incident  on  any  flame  is 
modeled  as  being  fully  absorbed  with  no  appreciable  effect  on  the  flame's 
radiating  characteristics. 

In  the  model  each  flame  is  represented  as  an  ensemble  of  point 
sources  above  the  nozzle,  along  its  axis  of  symmetry.  This  is  adequate  as 
long  as  enough  points  are  included  on  the  axis  to  allow  an  accurate 
integration,  and  if  the  effect  from  flame-flame  obscuration  is  properly 
estimated. 

The  parameters  of  the  inverted  cone  model  are:  (I)  flame  strength  or 
power,  measured  in  calories  per  second,  (2)  flame  width  or  diameter  (at  the 
"base"  of  the  inverted  cone),  measured  in  centimeters,  (3)  flame  height  at  the 
point  where  significant  radiation  seems  to  end,  also  in  centimeters,  (4)  and 
(5)  two  adjustable  parameters  that  quantitatively  represent  the  height 


variation  of  power  per  unit  length  divided  by  the  radius  of  the  flame  at  that 
height,  and  (6)  the  ground's  albedo.  These  six  parameters  are  discussed  more 
fully  in  Section  4. 


2-4  CALIBRATING  ALFGEE  MODELS. 

In  practice,  before  any  ALFGEE  model  is  used,  it  is  calibrated  to  the 
data.  A  calibration  is  defined  as  the  set  of  six  parameters  that  best 
represent  the  data  being  modeled,  where  this  set  will  probably  be  different 
from  the  set  for  another  physical  representation  of  the  flame's  geometry.  The 
right  cylinder  representation  was  the  one  most  widely  used  during  this  effort. 
It  was  also  the  one  used  to  predict  the  environments  for  MILL  RACE.  The  more 
recent  inverted  cone  model  had  not  been  exercised  much  during  this  effort. 

The  process  of  model  calibration  generally  proceeds  as  follows.  A 
set  of  parameters  are  chosen,  usually  based  on  past  experience.  Calculations 
for  points  where  data  exist  are  performed  and  these  compared  with  the 
experimentally  measured  values.  The  difference  between  the  predicted  and 
experimental  values  is  also  computed  and  is  called  the  residual.  The  absolute 
strength  parameter  of  the  flame  model  is  adjusted  to  give  a  zero  residual  sum. 
The  model's  measure  of  goodness  is  defined  to  be  the  smallness  of  the  sum  of 
the  variations,  where  a  variation  is  computed  as  the  residual  squared. 


The  entire  TRS  data  base  should  not  necessarily  be  used  when 

calibrating  a  model.  A  subset  should  be  chosen,  for  example,  excluding 
measurements  made  with  defective  calorimeters,  or  measurements  where  the 
aluminum  and  LOX  flow  rates  are  largely  different  from  current  practice.  The 

resultant  calibrations  can  be  compared  to  see  how  well  they  represent  the  data 

base.  It  is  important  to  inspect  individual  comparisons  to  be  sure 

differences  are  random  rather  than  biased. 


Since  ALFGEE  is  a  six  parameter  model,  any  subset  must  contain  at 
least  six  perfect  data  points  to  determine  a  unique  model.  Since  the  data 
taken  in  the  field  are  subject  to  variability,  much  more  data  must  be 


available  to  provide  a  statistically  meaningful  sampling. 


The  best  model  is  not  necessarily  that  with  a  zero  residual  sum  and 
the  lowest  sum  of  variations,  even  though  it  has  been  our  practice  to  select 
our  parameters  on  that  basis.  For  example,  if  a  comparison  of  the  data  with  a 
calibrated  model  showed  consistently  lower  flux  close  in  to  the  flames  and 
consistently  higher  flux  farther  out,  another  calibrated  model  might  be 
preferable  even  if  its  residual  sum  was  nonzero  and  it  had  a  larger  sum  of 
variations.  Forcing  the  residual  sum  to  be  zero  and  using  the  sum  of 
variations  as  the  measure  of  a  model's  ability  to  represent  some  data  base 
subset,  is  an  appropriate  technique  when  the  data  show  no  biases. 


2-5  EXAMPLES  OF  CALIBRATIONS. 

This  subsection  presents  results  of  four  different  calibration 
studies,  out  of  many,  that  have  been  performed  during  this  effort.  These  are 

examples  of  results  generated  with  plausible  model  parameters.  They  were  not 
produced  with  the  intent  to  derive  the  optimum  calibration.  Instead,  when 

someone  needed  a  prediction  within  a  short  time  constraint  and  only  a  meager 

set  of  data  existed,  we  would  run  the  necessary  calculations  to  quickly  arrive 
at  a  better  calibration,  and  then  provide  revised  model  estimates. 

Table  I  shows  comparisons  for  different  calibrations  of  the  right 

circular  cylinder  model.  The  table  consists  of  seven  columns  of  data.  The 
first  column  is  the  power  of  the  TRS  unit  (ABS.STR)  in  calories  per  second. 
The  second  (F.DIAM)  is  the  flame  diameter  in  centimeters.  The  third  (F.HGT) 

is  the  height  of  the  flame's  radiating  portion,  also  in  centimeters.  The 
fourth  (BK-HGT)  and  fifth  (BK-STR)  are  coefficients,  described  later,  that 
affect  the  variation  of  the  power  per  unit  length  (height)  of  flame  with  its 

height.  Both  coefficients  are  non-dimensional.  The  sixth  column  is  the 
ground  albedo,  another  model  parameter.  Finally,  the  seventh  column  is  the 
sum  of  the  variation  that  serves  as  our  measure  of  a  model's  goodness  of  fit. 

The  lower  this  number,  the  better  the  agreement. 


Table  1.  Variation  study  for  Navy  TRS  (no  passive  calorimeters). 


VARIATION  STUDY  FOR  NAVY  TRS  (NO  PASSIVE  CALORIMETERS) 


19-APR-82  IS  LAST  UPDATE 


DNA  TRS 


B.  CHAMBERS 


ABS.STR  F.DIAM  F.HGT  BK-HGT  BK-STR  ALBEDO  SUM-VAR 


For  the  chosen  subset,  representative  of  the  Navy  MILL  RACE  TRS  unit, 
large  variations  in  some  of  the  model  parameters,  as  illustrated  in  Table  I, 
did  not  make  much  difference.  Specifically,  it  shows  that  the  sum  of  the 
variances  is  insensitive  to  the  details  of  the  relationship  describing  how  the 
power  from  the  flame  depends  on  normalized  height.  The  two  parameters,  BK-HGT 
and  BK-STR,  quantify  this  relationship.  The  interpretation  that  best  fits 
this  case  is  that  the  models  are  equally  inadequate. 

In  the  above  example,  the  chosen  data  base  consisted  of  nineteen 
measurements  made  with  the  calorimeters  from  pre-MILL  RACE  test  firings  and 
the  active  ones  from  MILL  RACE.  Nine  of  the  measurements  were  taken  at  the 
test  firings  and  were  read  by  Miller  of  SAI.  The  nozzles  during  that  firing 
were  1.75  meters  above  ground.  Ten  of  the  measurements  were  made  on  MILL 
RACE.  Four  of  these  were  read  by  the  author.  Six  were  calculated  by  Bousek 
of  FCDNA  from  the  fluence  divided  by  2.4,  the  burn  time.  The  sum  of  the  flux 
measurements  for  the  nineteen  points  was  297  calories  per  square  centimeter 

per  second. 

Table  2  shows  the  results  for  the  same  Navy  TRS  unit  and  the  same 
nineteen  point  data  base  subset,  but  using  the  inverted  cone  model.  The 
columns  are  the  same  as  in  Table  I,  except  for  two  parameters  that  describe 
the  power  per  unit  flame  length  with  height  relatonship.  The  two  new 
parameters  are  the  luminosity  limit,  a  height  normalized  by  the  flame  height 

(LUM.LIM),  and  the  power  at  the  nozzle  (P@NOZZ),  which  is  relative  to  the 
power  per  unit  height  at  the  flame  height  (defined  as  unity).  This  latter 

parameter  is  normalized  to  conserve  energy  within  the  model. 

In  Table  2  the  sum  of  the  variation  is  much  more  sensitive  to 
variations  in  the  model  parameters.  The  inverted  cone  model  reduces  the  sum 
of  variances  by  fourteen  percent,  and  therefore  could  be  considered  a  better 
model  of  these  data. 

Table  3  compares  various  calibrations  of  the  inverted  cone  model  to 
another  subset  of  the  data  base.  Here  the  effect  of  adding  the  passive 

calorimeter  data  was  considered.  Twenty  two  additional  measurements  were 
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Table  2.  Variation  study  for  Navy  TRS  (does  not  include  passive  calors). 


VARIATION  STUDY  FOR  NAVY  TRS  (DOES  NOT  INCLUDE  PASSIVE  CALORS) 
21-APR-8 2  IS  LAST  UPDATE  DNA  TRS  B.  CHAMBERS 


ABS.STR 

F.DIAM 

F.HGT 

LUM.LIM 

P@NOZZ 

ALBEDO 

SUM-VAR 

30427800 

480.06 

640.08 

1.00 

1.00 

0.00 

129 

37174700 

480.06 

640.08 

1.00 

0.00 

0.00 

243 

32798800 

480.06 

640.08 

1.00 

0.50 

0.00 

160 

30427800 

480.06 

640.08 

1.00 

1.00 

0.00 

129 

27920100 

480.06 

640.08 

1.00 

2.00 

0.00 

105 

25807700 

480.06 

640.08 

1.00 

4.00 

0.00 

93 

24003900 

480.06 

640.08 

1.00 

10.00 

0.00 

89 

23111200 

480.06 

640.08 

1.00 

25.00 

0.00 

89 

22618100 

480.06 

640.08 

1.00 

100.00 

0.00 

90 

30427800 

480.06 

640.08 

1.00 

1.00 

0.00 

129 

28408500 

480.06 

640.08 

1.00 

1.00 

0.20 

198 

26640800 

480.06 

640.08 

1.00 

1.00 

0.40 

275 

25080200 

480.06 

640.08 

1.00 

1.00 

0.60 

356 

23692300 

480.06 

640.08 

1.00 

1.00 

0.80 

439 

22450000 

480.06 

640.08 

1.00 

1.00 

1.00 

521 

30427800 

480.06 

640.08 

1.00 

1.00 

0.00 

129 

34731200 

480.06 

640.08 

1.20 

1.00 

0.00 

155 

38761700 

480.06 

640.08 

1.40 

1.00 

0.00 

182 

50809200 

480.06 

640.08 

2.00 

1.00 

0.00 

255 

91850900 

480.06 

640.08 

4.00 

1.00 

0.00 

382 

29115100 

400.00 

640.08 

1.00 

1.00 

0.00 

135 

30427800 

480.06 

640.08 

1.00 

1.00 

0.00 

129 

33724300 

600.00 

640.08 

1.00 

1.00 

0.00 

147 

28636400 

480.06 

600.00 

1.00 

1.00 

0.00 

112 

30427800 

480.06 

640.08 

1.00 

1.00 

0.00 

129 

38487500 

480.06 

800.00 

1.00 

1.00 

0.00 

219 

18165500 

480.06 

400.00 

1.00 

25.00 

0.00 

153 

19939300 

480.06 

500.00 

1.00 

25.00 

0.00 

112 

22137400 

480.06 

600.00 

1.00 

25.00 

0.00 

92 

23111200 

480.06 

640.08 

1.00 

25.00 

0.00 

89 

24664700 

480.06 

700.00 

1.00 

25.00 

0.00 

90 

27452000 

480.06 

800.00 

1.00 

25.00 

0.00 

102 
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Table  3.  Variation  study  for  Navy  TRS  (includes  passive  calorimeters)  . 


VARIATION  STUDY  FOR  NAVY  TRS  (INCLUDES  PASSIVE  CALORIMETERS) 
18-APR-82  IS  LAST  UPDATE  DNA  TRS  B.  CHAMBERS 


ABS.STR 

FoDIAM 

F.HGT 

LUM.LIM 

P@NOZZ 

ALBEDO 

SUM- VAR 

28073100 

480.06 

640.08 

1.00 

1.00 

0.00 

1215 

33556300 

480.06 

640.08 

1.00 

0.00 

0.00 

1870 

30027600 

480.06 

640.08 

1.00 

0.50 

0.00 

1409 

28073100 

480.06 

640.08 

1.00 

1.00 

0.00 

1215 

25972000 

480.06 

640.08 

1.00 

2.00 

0.00 

1056 

24176300 

480.06 

640.08 

1.00 

4.00 

0.00 

960 

22622000 

480.06 

640.08 

1.00 

10.00 

0.00 

907 

21845900 

480.06 

640.08 

1.00 

25,00 

0.00 

891 

21415100 

480.06 

640.08 

1.00 

100.00 

0.00 

885 

28073100 

480.06 

640.08 

1.00 

1.00 

0.00 

1215 

26001800 

480.06 

640.08 

1.00 

1.00 

0.20 

1165 

24215200 

480.06 

640.08 

1.00 

1.00 

0.40 

1159 

22658300 

480.06 

640.08 

1.00 

1.00 

0.60 

1183 

21289500 

480.06 

640.08 

1.00 

1.00 

0.80 

1226 

20076700 

480.06 

640.08 

1.00 

1.00 

1.00 

1281 

28073100 

480.06 

640.08 

1.00 

1.00 

0.00 

1215 

31733800 

480.06 

640.08 

1.20 

1.00 

0.00 

1323 

35081500 

480.06 

640.08 

1.40 

1.00 

0.00 

1370 

45327400 

480.06 

640.08 

2.00 

1.00 

0.00 

1546 

81808900 

480.06 

640.08 

4.00 

1.00 

0.00 

1792 

26763100 

400.00 

640.08 

1.00 

1.00 

0.00 

1254 

39983700 

800.00 

640.08 

1.00 

1.00 

0.00 

916 

23188900 

480.06 

500.00 

1.00 

1.00 

0.00 

972 

34628200 

480.06 

800.00 

1.00 

1.00 

0.00 

1552 
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added  and  then  some  of  the  calibrations  in  the  previous  table  were  redone.  The 
variability  of  these  additional  data,  as  measured  by  the  sum  of  the  variations 
divided  by  the  number  of  points,  increased.  The  variability  is  increased  by 
adding  the  passive  calorimeter  data.  The  added  data  were  calculated  from  the 
measured  fluence  by  Bousek  of  FCDNA.  The  approach  of  using  measured  fluences 
to  determine  flux  can  be  misleading  in  the  opinion  of  this  writer. 

Table  4  presents  calibrations  to  data  taken  from  two  sets  of  nozzles 
using  the  inverted  cone  model.  The  two  sets  were  lumped  together  because  they 
were  both  spaced  seven  feet  apart.  The  two  sets  were  the  UK- 1  and  UK-2  MILL 
RACE  nozzles.  UK-1  was  an  eight  nozzle  array  and  UK-2  was  a  four  nozzle 
array.  The  data  consisted  of  thirty  three  measurements:  ten  were  the  UK- 1 
pre-MILL  RACE  test  firings  for  UK-1  and  were  read  by  Miller  of  SAI;  eleven 
were  from  the  pre-MILL  RACE  test  firings  for  UK-2  and  were  read  by  Lottery  of 
SAI;  eight  were  the  UK- 1  MILL  RACE  results  read  by  this  writer;  and  four  were 
the  UK-2  MILL  RACE  results  also  read  by  this  writer.  The  sum  of  the  measured 
values  was  794  calories  per  square  centimeter  per  second.  The  lowest  sum  of 
variances  was  323  in  flux  squared  units.  This  result  was  not  much  different 
from  the  one  obtained  by  using  the  best  modeling  for  the  Navy  unit  (Table  2), 
the  latter  shown  at  the  bottom  of  Table  4  for  the  UK  data.  The  residual  is  29 
calories  per  square  centimeter  per  second,  and  is  non-zero  because  the 
parameters  were  deliberately  not  changed.  This  latter  comparison  led  us  to 
believe  that  it  is  possible  to  build  a  unified  model,  that  is,  one  with 

parameters  that  are  independent  of  nozzle  spacing. 

2-6  DNA  TRS  DATA  BASE. 

A  data  base  has  been  compiled  from  the  experimental  measurements 

recorded  during  the  recent  years  of  TRS  firings.  Any  subset  of  the  data  base 
can  potentially  be  used  to  calibrate  any  ALFGEE  model,  e.g.  cylindrical  or 
inverted  cone.  Depending  upon  the  success  achieved  in  replicating  the 
experimental  data,  appropriate  subsets  of  the  data  base  can  provide  an 
indication  of  the  overall  variability  of  the  data.  Whether  any  variation  hos 

been  caused  by  experimental  error  or  instead  by  lack  of  repeatability  is  not 
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Table  4.  Variation  study  for  UK's  TRS  (does  not  include  passive  calors)  . 


VARIATION  STUDY  FOR  UK'S  TRS  (DOES  NOT  INCLUDE  PASSIVE  CALORS) 
26-APR-8 2  IS  LAST  UPDATE  DNA  TRS  B.  CHAMBERS 


ABS.STR 

F.DIAM 

F.HGT 

LUM.LIM 

P@NOZZ 

ALBEDO 

SUM-VAR 

24311200 

480.06 

640.08 

1.00 

25.00 

0.00 

371 

20369200 

480.06 

350.00 

1.00 

25.00 

0.00 

858 

20596000 

480.06 

400.00 

1.00 

25.00 

0.00 

670 

21073600 

480.06 

450.00 

1.00 

25.00 

0.00 

536 

21724500 

480.06 

500.00 

1.00 

25.00 

0.00 

447 

23445900 

480.06 

600.00 

1.00 

25.00 

0.00 

379 

25754100 

480.06 

700.00 

1.00 

25.00 

0.00 

368 

28395300 

480.06 

800.00 

1.00 

25.00 

0.00 

395 

33920700 

480.06 

640.08 

1.00 

1.00 

0.00 

620 

25426300 

480.06 

640.08 

1.00 

10.00 

0.00 

381 

24311200 

480.06 

640.08 

1.00 

25.00 

0.00 

371 

23701300 

480.06 

640.08 

1.00 

100.00 

0.00 

368 

24311200 

480.06 

640.08 

1.00 

25.00 

0.00 

371 

23069900 

480.06 

640.08 

1.00 

25.00 

0.20 

424 

2194930 0 

480.06 

140.08 

1.00 

25.00 

0.40 

491 

21370600 

300.00 

640.08 

1.00 

25.00 

0.00 

334 

22019700 

350.00 

640.08 

1.00 

25.00 

0.00 

342 

22835100 

400.00 

640.08 

1.00 

25.00 

0.00 

339 

23728800 

450.00 

640.08 

1.00 

25.00 

0.00 

360 

24311200 

480.06 

640.08 

1.00 

25.00 

0.00 

371 

27017000 

600.00 

640.08 

1.00 

25.00 

0.00 

374 

24311200 

480.06 

640.08 

1.00 

25.00 

0.00 

371 

25351200 

480.06 

640.08 

1.50 

25.00 

0.00 

354 

26441300 

480.06 

640.08 

2.00 

25.00 

0.00 

342 

29064100 

480.06 

640.08 

3.00 

25.00 

0.00 

348 

31615000 

480.06 

640.08 

4.00 

25.00 

0.00 

340 

24828600 

400.00 

640.08 

2.00 

25.00 

0.00 

323 

23111200 

480.06 

640.08 

1.00 

25.00 

0.00 

391 

*  NOTE:  THIS  IS  A  SPECIAL  CASE  (10.5)  SUM  OF  RESIDUAL>  29 


easily  deduced. 


Subsets  of  the  data  base  have  been  chosen  (and  used  to  establish 
different  calibrations)  as  statistically  meaningful  collections  representing 
either:  (I)  a  particular  type  of  event,  such  as  the  collection  of  all  multi¬ 

nozzle  events  wherein  the  nozzles  are  separated  by  7  feet,  or  (2)  where 
certain  suspect  measurements  have  been  removed  for  example  for  the  3.5  feet 
nozzle  separation  data.  The  most  suspicious  data,  those  obtained  from  early 
experiments  performed  with  different  reactant  flow  rates,  were  removed  in  an 
attempt  to  avoid  skewing  the  data. 

Appendix  B  presents  some  of  the  peak  flux  data  taken  during  the 
development  program,  as  reported  to  this  author,  along  with  comparisons  of  an 
early  model  based  on  the  right  cylinder.  The  recorded  peak  value  has  been 

found  to  depend  on  who  is  reading  the  oscilliscope  traces,  primarily  due  to 

the  interpretation  of  what  is  peak  flux.  Some  take  the  absolute  peak,  they 
don't  average  over  the  rapid  oscillations.  Others  use  a  faired  curve  through 
the  data  removing  the  high  frequencies  (kiloHertz  and  above).  Therefore,  the 
appendix  includes  more  comparisons  than  would  otherwise  be  needed. 

2-7  CALORIMETER  LOCATIONS. 

Figures  I  through  6  illustrate  where  the  calorimeters  have  been 
located  relative  to  the  flames.  Each  calorimeter  is  represented  by  a  small 
circle,  which  is  connected  to  its  projection  on  the  ground.  The  circle  is 
drawn  slightly  to  the  calorimeter's  right.  The  flames  are  pictorally  shown  as 
four  straight  lines  always  starting  at  the  nozzle.  Each  of  these  figures  are 
projections;  they  are  perspective  only  for  an  eye  position  infinitely  far 
away.  The  coordinate  scale  is  shown  in  the  upper  left  hand  corner  and  is  in 
centimeters.  This  scale  refers  to  the  distance  at  the  end  of  each  coordinate 

axis.  Each  axis  is  subdivided  into  ten  intervals,  each  interval  bounded  by 

tic  marks. 

Note  that  the  pre  MILL  RACE  test  firings  of  both  the  UK-1  eight 
nozzle  array  and  the  Navy  four  nozzle  array  were  positioned  above  the  ground 


since  both  arrays  were  too  long  to  fit  in  the  TRS  pit  at  the  KAFB  site.  The 
nozzles  were  above  the  y-axis  in  the  figures  and  were  at  a  height  of  1.75 
meters.  Referring  to  Figure  4,  this  height  of  1.75  meters  when  divided  by  the 
coordinate  scale,  7.54  meters,  produces  a  value  of  0.232,  which  corresponds  to 
2.32  tic  marks  up  the  z  axis.  This  illustrates  how  to  use  the  coordinate 
scale. 

2-8  BIAS  OF  DATA  BASE  WITH  EVENT. 

One  of  the  better  calibrations  was  chosen  to  be  compared  with  all  the 
data  available  in  the  data  base  (including  MILL  RACE  diagnostic  calorimeters) 
to  find  any  evidence  of  bias.  The  parameters  for  the  chosen  model  are  shown 
above  in  Table  4  as  the  last  entry.  That  calibration  was  to  the  Navy  data, 
which  did  not  include  the  passive  calorimeters. 

The  ratio  of  each  experimental  datum  over  the  model  value  was 
computed  and  stored  in  a  list  along  with  its  x  (range)  and  z  (height) 
coordinates.  Each  calorimeter  and  event  identification  was  also  included. 
Appendix  C  contains  two  tables  that  present  this  information,  which  is  based 
on  the  data  available  to  this  writer  up  to  and  including  MILL  RACE,  in  two 

different  orders.  The  first  table  shows  the  list  in  chronological  order.  The 
second  table  shows  this  list  sorted  on  the  value  of  the  result  ratio.  It  is 

included  as  Appendix  C  for  reference. 

Figures  7  and  8  show  the  variation  of  the  ratio  of  measured  values  to 
calculated  ones  for  two  sets  of  data.  The  first  set  is  for  the  TRS  unit  at 
the  DNA  site  at  KAFB.  Initially  this  consisted  of  two  nozzles,  or  one  half  of 
a  normal  unit,  but  later  was  modified  to  include  two  more  nozzles.  The  second 
is  for  the  four  sets  of  TRS  units  built  specifically  for  use  on  MILL  RACE. 

There  does  not  seem  to  be  any  discernible  bias  with  event  nor  with  time. 
However,  there  is  clearly  considerable  variation  in  how  well  the  data  agree. 

The  distribution  of  variability  with  ratio  is  shown  as  Figure  9.  In 

general  since  the  ratio  is  used,  figures  of  this  type  are  more  useful  when  the 
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Figure  9.  Distribution  of  ratio  differences. 


log  of  the  ratio  is  shown.  The  information  to  be  gained,  however,  did  not  seem 
significant  enough  to  replot  the  data.  What  this  figure  does  show  is  that  in 
general  the  spread  in  data  appears  almost  normally  distributed. 

Taking  the  data  shown  in  Figures  7  and  8,  subtracting  one  from  the 
ratio  and  then  accumulating  the  result  as  datum  number  increases,  Figures  10 
and  I  I  are  generated.  If  there  was  no  bias  with  time,  these  figures  would 
oscillate  about  zero.  Where  there  is  a  positive  or  negative  slope  over  a 
number  of  data  points,  an  event  is  indicated  which  on  average  is  higher  or 
lower,  respectively,  than  the  model  calculations.  One  such  dominant  slope 
region  is  on  Figure  10  from  datum  number  45  to  66.  This  region  corresponds 
exactly  to  Events  D-30,  D-31,  D-32,  D-33,  D-34,  and  D-35.  Whether  this  same 
region  may  have  been  seen  in  Figure  7  is  debatable.  The  region  is  quite 
discernible  in  Figure  10,  which  makes  the  latter  type  of  presentation  useful 
for  data  analysis.  The  trend  is  also  seen  in  Figure  II  for  the  UK- 1  and  BRL 
MILL  RACE  pretest  firings. 

2-9  BIAS  OF  DATA  BASE  WITH  RANGE  AND  HEIGHT. 

In  order  to  establish  whether  there  existed  a  bias  with  calorimeter 
field  position,  the  figures  within  this  subsection  were  prepared.  All  the 
active  calorimeter  data  were  separated  into  two  groups.  One  group  was  for  the 
DNA  TRS  site  unit,  and  the  other  group  were  all  the  MILL  RACE  units.  Three 
figures  were  prepared  for  each  set.  The  first  shows  the  projection  of  the 
locations  of  each  calorimeter  on  the  range-height  plane.  The  second  shows  the 
ratio  of  the  measured  value  to  the  computed  one  for  the  calibration  used  in 
the  preceding  subsections  as  a  function  of  range.  The  third  shows  this  same 
ratio  as  a  function  of  height. 

For  the  DNA  TRS  site  unit,  all  data  from  Events  D-15  to  Pentolite  5 
are  shown  in  figures  12  through  14.  The  first  thing  that  should  be  noticed 
is  that  the  field  has  only  been  modestly  covered.  However,  the  nozzles  have 
been  near  the  ground  in  all  cases.  The  second  thing  to  be  noticed  is  that 
there  does  not  appear  to  be  a  clear  bias  of  the  variability  with  range. 
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However,  there  does  appear  to  be  a  slight  bias  as  a  function  of  calorimeter 
height.  This  bias  with  height  is  looked  at  more  closely  in  the  following 
subsection. 

Figures  15  through  17  show  similar  data,  but  instead  for  the  MILL 
RACE  nozzles.  Here  again  there  does  not  appear  to  be  any  overall  bias, 
neither  with  range  nor  with  height.  The  interpretation  in  this  case  is  not  as 
straightforward,  because  during  two  of  the  pretest  firings,  the  nozzles  were 
1.75  meters  above  the  ground. 


2-10  A  CLOSER  LOOK  AT  BIAS  OF  DATA  BASE  WITH  HEIGHT. 

Some  of  the  data  presented  in  the  previous  subsection  showed  an 

apparent  bias  of  the  ratio  as  a  function  of  height.  In  this  subsection, 

smaller  samples  of  the  data  base  are  presented  that  were  analyzed  to  see  if 

the  height  bias  can  be  made  more  visible,  or  if  it  is  a  function  of  event,  or 
equivalently,  time. 

Figure  18  shows  a  plot  of  ratio  versus  height  for  all  the  data 
earlier  than  the  Pentolite  series,  that  is,  Events  D-15  to  D-45.  A  slightly 
better  hint  of  a  height  bias  is  seen.  The  alternative  plot  style  also  makes 

it  clear  where  the  points  are  really  located.  Subdividing  these  data  further 
into  two  subsets,  yields  Figures  19  and  20.  The  choice  of  Event  D-26  as  the 
break  point  was  somewhat  arbitrary.  Nevertheless,  we  see  that  the  height  bias 
is  very  clear  for  the  very  early  data,  whereas  it  is  not  clear  that  a  bias 
exists  for  the  later  data. 

During  the  early  testing,  considerable  work  was  being  done  to  improve 
the  flame  output.  In  fact,  the  reason  the  data  earlier  than  D-15  has  not  been 
used  in  the  modeling  effort  is  that  reasonable  operating  characteristics  had 
not  been  established. 
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A  CLOSER  LOOK  AT  REPRODUCIBILITY. 


2-1 1 

After  having  inspected  the  entire  data  base  taken  up  through  MILL 
RACE,  having  gained  a  certain  confidence  in  the  data,  and  having  finished  the 
time  urgent  work  for  MILL  RACE,  we  went  back  to  look  at  the  flux  data  to  see 
if  we  could  arrive  at  some  measure  of  reproducibility.  The  most  readily 
comparable  data  to  establish  this,  as  it  turns  out,  are  the  early  data.  They 

are  unique  in  that  the  calorimeters  were  never  moved,  except  for  FX-I  after  it 
was  established  that  both  sides  of  the  TRS  unit  generated  nearly  identical 
flux  levels  at  the  same  range.  Furthermore,  the  TRS  unit  consisted  of  the 
same  two  nozzles. 

Figures  21  through  24  show  the  peak  flux  measurements  for  the  four 
calorimeters  that  were  available  in  the  early  test  program.  The  calorimeters 

were  labeled  as  FX-I,  FX-2,  FX-3  and  FX-4.  The  position  for  calorimeter  FX-I 
changes  immediately  after  event  D-21.  Clearly,  the  data  recorded  for  event 
D-27  is  anomalous. 

Replotting  these  data,  leaving  out  event  D-27  and  splitting  the  plot 

for  FX-I  to  reflect  the  two  positions,  more  clearly  reflects  the  true  nature 

of  the  data.  These  are  shown  as  Figures  25  through  29.  The  mean  for  FX-I  for 
events  D-15  through  D-21  is  41.3  calories  per  square  centimeter  per  second. 
The  corresponding  unbiased  standard  deviation,  s,  is  1.39  in  flux  units.  This 
compares  quite  well  with  the  mean  for  FX-4  for  events  D-15  through  D-31  with 
D-27  removed.  The  mean  for  FX-4  is  40.9  and  its  unbiased  standard  deviation 
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Statistics  for  Events  D-15  through  D-31 


Calorimeter 

Height 

(ft) 

Event 

Start  End 

Mean 

m 

Standard 

Deviation 

s/m 

FX-! 

5.0 

D-15 

D-21 

41.3 

1.39 

FX-I 

7.5 

D-23 

D-31 

59.9 

3.84 

FX-2 

10.0 

D-15 

D-31 

13.5 

2.37 

0.18 

FX-3 

10.0 

D-15 

D-31 

31.9 

5.44 

FX-4 

5.0 

D-15 

D-31 

40.9 

2.90 

0.07 

The 

flux  varies  the 

most  at 

the  higher  calorimeters.  This 

is  where 

the  effects 

of  wind 

will  be 

most  pronounced. 

From  these 

data  we  get  an 

indication  that  the 

shot  to 

shot  repeatability  will  depend 

on  where  the 

measurement 

is  being 

made.  Within  the 

limited  volume  of  space  where 

the  data 

were  taken, 

the  data 

suggest 

that  the 

standard 

deviation  is  about  7 

percent 

below  5  feet 

(1.52  meters)  in 

height  for  ranges 

within  6  feet 

(1.83  meters)  of 

the  nozzles, 

and  is  18  percent 

at  10  feet  (3.05 

meters)  height 

within 

a  range 

of  6  feet  (1.83  meters)  to  12  feet  (3.66  meters). 


SECTION  3 


COMPARISONS  WITH  MILL  RACE 

This  section  presents  comparisons  of  various  ALFGEE  models  with 
calorimeter  data  taken  on  MILL  RACE.  These  comparisons  illustrate  the  size  of 
the  differences  that  can  be  expected  as  a  result  of  calibrating  ALFGEE  to 
limited  data.  Twenty  diagnostic  calorimeters  were  fielded  by  Ballistics 
Research  Laboratory  (BRL)  for  DNA.  One  of  the  calorimeters  failed  to  produce 
any  data.  About  half  of  the  rest  agreed  very  well  with  the  model  predictions. 
Two  calorimeters  gave  very  different  values  from  what  was  expected.  One  of 
them  was  about  50  percent  too  high  and  the  other  was  about  50  percent  too  low. 
On  the  whole,  the  modeling  provided  a  reasonably  accurate  prediction 
capability. 

Not  all  of  the  TRS  units  burned  as  had  been  hoped.  One  in 
particular,  UK- 1,  had  two  nozzles  that  did  not  ignite,  because  a  jammed  flow 
valve  stopped  aluminum  flow.  Furthermore,  another  nozzle  on  this  same  unit 
apparently  did  not  burn  well  due  to  a  faulty  pilot  flame.  Based  on 
performance  of  the  TRS  units  on  MILL  RACE,  experimenters  should  consider 
possible  deviations  from  planned  environments,  especially  at  locations  where 
little  experimental  data  exist  upon  which  to  base  the  modeling. 


3-1  ALFGEE  MILL  RACE  COMPARISONS  (CYLINDRICAL). 

The  comparisons  of  the  peak  flux  data  taken  on  Ihe  BRL/HDL  HATS  TRS 
is  summarized  first  since  its  nozzle  separation  was  the  same  as  all  the  data 
taken  prior  to  the  pretest  firings  of  the  MILL  RACE  TRS  units.  Its  nozzles 
were  spaced  3.5  feet  apart.  The  agreement  with  the  model,  as  it  existed  prior 
to  MILL  RACE,  is  shown  in  Table  5  (Model  A  cylindrical).  Only  three  out  of 
four  diagnostic  calorimeters  yielded  data  for  this  four  nozzle  burn.  Only  two 


Table  5.  Comparisons  of  BRL  TRS  at  MILL  RACE  with  DNA  ALFGEE  model 


Peak  Flux  (cal/cm/cm/s) 


Mill  Race 

Model-A 

Model-B 

Model-82 

F-l 

62.0 

39.8 

39.6 

41.0 

F-2 

50.6 

55.6 

52.2 

54.8 

F-3 

No  data 

55.9 

52.5 

55.1 

F-4 

37.6 

40.7 

40.4 

41.8 

Note:  Model-A  is  best  model  used  prior  to  MILL  RACE. 
Model-B  is  recalibration  to  162  point  data  base. 
Model-X  is  inverted  cone  3.5'  calibration. 


of  these  calorimeter  records  were  consistent  with  the  expected  environments. 
The  cylindrical  model  was  calibrated  to  the  meager  data  obtained  during  the 

BRL  TRS  unit  pretest  firings.  This  calibration  improves  the  agreement  of  model 
and  MILL  RACE  data.  We  assert  that  these  comparisons  show  that  the  data  from 
calorimeter  F-l  is  highly  suspicious. 

The  results  for  UK-2  are  shown  in  Table  6,  where  the  calibration 

based  on  the  sparse  eleven  point  7  feet  spacing  data  base  yields  the  best 

agreement.  Calorimeter  F-l  is  the  only  one  that  disagrees  by  more  than  three 
per  cent  from  the  expected  environment. 

The  results  for  UK- 1,  the  only  8  nozzle  system  (7'  spacing),  was 

unique  in  that  two  nozzles,  numbers  3  and  4,  did  not  inject  any  aluminum  into 

the  air.  Not  only  did  these  nozzles  fail  to  contribute  any  flux  to  the 

environment,  but  the  fact  that  they  did  not  ignite  makes  it  difficult  to 
select  an  appropriate  calibration  with  which  to  measure  the  performance  of  the 
TRS  system.  Furthermore,  the  films  of  the  BRL  TRS  shown  during  the  MILL  RACE 
symposium  show  that  nozzle  7  burned  in  an  unusual  manner  and  may  have 
contributed  significantly  less  flux  than  otherwise  would  have  been  expected. 

Because  the  effect  of  the  nozzles  either  not  firing  or  misfiring  was 

uncertain,  a  number  of  calculations  with  different  models  were  performed  with 
the  intent  to  estimate  the  magnitude  of  the  possible  effects.  The  results, 
summarized  in  Table  7,  show  that  the  calibration  to  the  1 1  point  7'  spacing 
data  again  yielded  the  best  overall  agreement  but  only  when  nozzles  3  and  4 
are  turned  off.  No  evaluation  of  the  effect  of  additionally  turning  off 
nozzle  7  has  been  made. 

Finally,  the  10.5'  spacing  Navy  4  nozzle  TRS  results  are  summarized 
in  Table  8.  A  study  was  performed  to  determine  whether  any  use  of  the  data 
obtained  during  the  experimental  calibration  of  the  Navy  TRS  could  be  made  for 
calibrating  a  model  based  on  a  cylindrical  flame.  It  was  possible  to 

substantially  improve  the  agreement  for  the  two  calorimeters  that  were  placed 
farthest  out.  However,  the  results  for  the  two  close-in  calorimeters  were 
degraded. 


Table  7.  Comparisons  of  UK- 1  TRS  at  MILL  RACE  with  DNA  ALFGEE  model. 


Table  8.  Comparisons  of  NAVY  TRS  at  MILL  RACE  with  DNA  ALFGEE  model. 


Peak  Flux  (cal/cm/cm/s) 


Mill  Race 

Model-A 

Model-C 

Model-D 

Model-Z 

Mode  1-82 

F-l 

13.2 

15.4 

14.3 

16.5 

15.2 

15.4 

F-2 

13.2 

9.7 

1  1.4 

13.2 

1  1.9 

10.0 

F-3 

12.0 

9.7 

II. 1 

12.8 

1 1.6 

10.0 

F-4 

12.0 

18.8 

17.4 

20.1 

18.5 

18.9 

Model-A 

is  model  used 

prior  to 

MILL  RACE 

Model-C  is  recalibration  to  I  I  point  7'  data  base. 
Model-D  is  recalibration  to  9  point  10.5'  data  base. 
Model-Z  is  Inverted  calibration  to  19  pt  10.5'  data. 


It  should  be  mentioned  that  the  calibration  shot  performed  on  the 
Navy  TRS  unfortunately  included  only  calorimeters  placed  at  the  same  height  of 


3.2  meters.  That  particular  arrangement  makes  it  very  difficult  to  calibrate 
any  model  with  any  degree  of  confidence. 

3-2  USE  OF  INVERTED  CONE  MODEL  WITH  TRS  DATA  BASE. 

Predictions  based  on  the  new  inverted  cone  model  with  the  changed 
view  angle  algorithms  were  compared  with  the  entire  TRS  data  base  to  determine 
whether  it  is  indeed  an  improvement,  and  what  are  the  best  parameters  to  model 
the  experimental  results.  The  conclusion  thus  far  seems  to  be  that  all  of  the 
recent  data  could  be  represented  with  one  calibration,  instead  of  having  to 
resort  to  a  different  calibration  for  each  nozzle  spacing.  In  other  words,  it 
may  be  possible  to  explain  the  effects  of  greater  nozzle  separation  on 
geometric  and  obscuration  considerations.  If  in  fact  this  conclusion  is 
finally  supported,  a  possibly  sound  basis  for  better  quantifying  the  shot  to 
shot  variability  may  now  exist. 

As  further  testimony  to  the  apparent  improvement  gained  with  the  new 
inverted  cone  model,  we  would  like  to  mention  that  two  of  the  model  parameters 
were  taken  directly  from  cursory  measurements  on  still  pictures  of  flames  at  a 
time  of  peak  flux.  This  latter  fact  greatly  improves  our  confidence  that  the 
newest  model  may  approximate  very  well  the  TRS  flame  thermal 
radiation  environment. 

3-3  COMPARING  CALIBRATIONS  FOR  MILL  RACE. 

A  number  of  comparisons  between  pairs  of  calibrations  were  performed 
to  evaluate  the  sensitivity  of  the  agreement  with  the  MILL  RACE  data  to  model 
parameters.  Some  of  the  comparisons  were  useful;  many  were  not.  In  order  to 
acquaint  the  reader  with  the  relative  changes  made  with  different  model 
parameters,  the  following  comparison  is  presented. 


Table  9  is  one  example  of  the  comparisons  made  between  different 
models.  This  table  shows  four  calibrations  and  compares  the  results  of  each 
model  prediction  for  the  set  of  data  that  includes  both  the  UK- 1  and  UK-2 
units  pre  MILL  RACE  test  firings  and  the  MILL  RACE  firings.  Since  two  of  the 
nozzles  did  not  fire  for  the  UK-1  unit  on  MILL  RACE,  their  source  strength  was 
set  to  zero  in  the  model  predictions.  No  account,  however,  was  taken  for  the 
badly  fired  nozzle,  since  we  had  no  way  of  knowing  what  its  strength  was  on 
the  test. 

Overall  the  comparisons  presented  here  show  that  the  cylindrical 
model  gave  slightly  better  results  than  the  three  calibrations  to  the  inverted 
cone  model.  While  this  will  sometimes  be  the  result  for  these  types  of 
comparisons,  the  general  result  has  been  that  the  inverted  cone  model  can  more 
consistently  agree  with  experimental  data.  This  particular  comparison 
illustrates  that  choosing  the  best  geometrical  representation  for  modeling  is 
not  necessarily  straightforward. 

3-4  COMPARING  FLUX  FIELD  MAPS. 

Substantial  differences  can  exist  between  the  predictions  of  one 
model  over  another  at  some  locations.  Less  often,  this  can  even  be  true  if 
both  were  calibrated  to  similar  data.  It  happens  quite  frequently,  however, 
when  the  models  have  to  be  used  in  regions  where  the  results  must  be 
extrapolated. 

The  following  tables,  10  through  13,  show  some  comparisons  between 
pairs  of  two  models  that  shed  some  light  on  where  biases  can  be  significant. 
For  example,  the  comparison  of  model  A  and  model  D  for  the  Navy  MILL  RACE  TRS 
unit  shows  similar  results  at  the  ground,  but  model  D  predicts  a  significantly 
higher  flux  (about  40  percent)  at  3  meters  height. 

These  types  of  comparisons  can  be  important  for  planning  HE  tests 
because  they  illustrate  the  uncertainties  in  model  predictions.  When  these 
uncertainties  are  large  and  when  this  is  important  to  the  experimenter,  he 


Table  9.  Comparisons  of  models  for  UK- 1  &  2. 


Table 


Table  10.  Navy  model  comparison  for  MILL  RACE  maps 


NAVY  MODEL  COMPARISON  FOR  MILL  RACE  MAPS  29DEC81 


Table  11.  BRL  model  comparisons  for  MILL  RACE  maps. 


BRL  MODEL  COMPARISONS  FOR  MILL  RACE  MAPS  29DEC81 


Table  12.  UK-2  model  comparisons  for  MILL  RACE  maps. 
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should  insist  on  more  calibration  of  the  TRS  unit  to  be  used. 


3-5  A  NOTE  ON  THE  MILL  RACE  DATA. 

Comparison  of  experimental  measurements  with  theoretical  predictions 
or  model  results  should  include  a  presentation  of  the  experimental  data.  Since 
the  data  are  still  unpublished,  a  preliminary  result  for  calorimeter  F-4 
exposed  to  the  BRL  unit  on  MILL  RACE  is  shown  as  Figure  30. 

Calorimeter  F-2  exhibited  a  waveform  similar  to  that  shown  in  Figure 
30  for  F-4.  Calorimeter  F-l,  on  the  other  hand  showed  a  significantly 
different  waveform.  Unfortunately,  calorimeter  F-3  failed. 

All  four  calorimeters  had  been  placed  at  the  same  range,  1.9  meters, 
from  the  nozzles  on  the  ground  zero  side  of  the  array.  They  also  had  been 
located  at  the  same  height  of  2.0  meters.  One  calo'  meter  was  placed  in  front 
of  each  of  the  four  nozzles  in  an  effort  to  record  each  nozzles’  performance. 
Since  the  field  of  view  of  each  calorimeter  was  not  limited  some  contribution 
from  the  other  nozzles  was  measured  at  each  calorimeter  (those  that  worked, 
that  is). 


The  agreement  of  the  BRL  MILL  RACE  data  with  the  model  calculations 
was  good  (within  ten  percent)  for  calorimeters  F-2  and  F-4.  Since  these  were 
the  ones  with  similar  waveforms,  the  agreement  is  encouraging.  The  fact  that 
the  agreement  is  poor  for  F-l  does  not  immediately  lead  to  a  conclusion.  Did 
wind  affect  the  results  for  this  calorimeter?  Is  the  model  reliable?  Until 
the  the  number  of  diagnostic  calorimeters  is  increased  on  these  tests,  many 
questions  will  remain  unanswered.  If  the  variability  of  the  experiment  was 
substantially  less,  then  it  would  be  sufficient  to  use  just  a  few  diagnostic 
calorimeters.  Since  that  is  not  the  case,  it  should  be  our  practice  to  use  at 
least  three  calorimeters  per  nozzle. 
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waveform:  BRL  calorimeter  F-4 


SECTION  4 


THE  ALFGEE  MODELS 

This  section  presents  the  ALFGEE  algorithms  used  in  the  modeling 

effort.  For  those  readers  interested  in  a  brief  introduction  to  the  pertinent 

parts  of  radiation  theory,  we  suggest  a  review  of  pages  183-188  in  Reference 

4.  In  addition,  pages  36  and  82  of  Reference  I  present  the  equations  for 

calculating  the  flux  from  either  an  extended  source  or  an  ensemble  of  point 
sources.  This  review  is  recommended  for  those  that  wish  to  follow  the  ensuing 
discussion. 


4-1  RADIATION  FROM  THE  FLAMES. 

The  calculation  of  the  radiant  flux  across  a  surface  facing  the 
flames  is  simply  an  integration  over  the  sources  of  radiation,  accounting  for 
the  effects  of  geometry.  A  surface  integral  is  evaluated  over  an  extended 
source.  Addition  of  all  the  contributions  from  each  source  is  performed  for 
an  ensemble  of  point  sources. 

An  extended  source  can  be  approximated  by  a  set  of  point  sources,  and 
then  the  flux  computed  as  for  an  ensemble  of  point  sources.  The  key  is 
choosing  where  to  place  the  points.  This  is  not  much  different  from 
numerically  solving  a  surface  integral.  In  fact,  to  obtain  the  same  accuracy 
in  both  methods  would  require  about  the  same  amount  of  calculation.  There  is 
essentially  no  difference  between  treating  the  surface  as  point  sources  and  as 
small  surfaces  when  the  extended  source  is  opaque. 

When  certain  symmetries  exist  it  is  possible  to  replace  a  surface 
integral  with  a  single  point.  An  example  of  this  is  radiation  from  an 
isotropic  radiating  sphere  out  to  a  collecting  surface,  which  is  outside  of 
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the  sphere  and  is  facing  the  sphere.  As  long  as  the  entire  sphere  that  can  be 
possibly  visible  to  a  collecting  surface  (depends  on  range)  is  visible  to  it, 
the  received  radiation  would  be  the  same  as  that  from  a  point  source  placed  at 
the  center  of  the  sphere. 

For  modeling  purposes,  we  have  chosen  to  approximate  the  extended 
source  as  a  set  of  point  sources  placed  on  the  axis  of  symmetry  of  the  flame. 
The  intent  is  to  avoid  integration  over  the  width  of  the  flame,  thus  keeping 
the  amount  of  calculation  at  a  minimum.  This  approach  introduces  small 
errors;  these  are  unimportant  because  the  model  is  ultimately  calibrated  to 
data.  Hence,  any  bias  introduced  with  height  will  tend  to  be  removed  with  the 
subsequent  calibration.  This  approach  is  an  outgrowth  of  the  earlier  work  on 
the  bag  source  reported  in  Reference  I.  Similarily,  the  criterion  of  how  many 
point  sources  need  to  be  distributed  along  this  line  was  taken  from  this 
earlier  work.  An  additional  constant  factor  was  included  to  ensure  accuracy 
in  unusual  geometries. 

Once  it  has  been  established  how  many  points  are  needed  to  represent 
a  flame  (this  is  done  for  each  field  point),  the  source  strength  of  the  entire 

flame  is  distributed  in  proportion  to  the  power-per-unit-length-model,  which 
is  described  later.  Since  the  number  of  points  needed  to  represent  each  flame 
changes  with  field  position,  the  source  strength  of  each  point  source  is 

normalized  to  the  total  from  the  ensemble. 

Features  such  as  normalization  help  ensure  a  continuity  in  the 
results  when  the  collecting  surface  moves  about  in  the  illuminated  field.  For 
example,  as  the  field  point  is  moved  toward  the  flames,  the  number  of  source 

points  increases  as  the  range  decreases,  yet  the  computed  flux  apparently 
exhibits  no  discontinuity.  Naturally,  these  smooth  transitions  should  not  be 
important  to  the  experimenter  nor  required  by  him.  Nevertheless,  they  do  help 
reduce  noise  in  the  computed  results,  which  is  useful  when  drawing  contours. 

The  point  sources  are  distributed  along  the  centerline  of  a  nozzle 
and  span  the  height  of  the  radiating  part  of  the  flame.  The  height  of  the 
flame  is  adjustable  by  the  analyst.  In  early  models,  the  resemblance  of  this 


parameter  to  an  actual  flame  height  was  coincidental.  More  recently,  however, 
the  best  fits  using  the  latest  inverted  cone  model  seem  to  yield  a  "flame 

height"  parameter  that  appears  to  agree  well  with  photographs. 

4-2  ALF GEE  COORDINATE  SYSTEM. 

Figure  31  presents  the  ALFGEE  coordinate  system  used  throughout  this 
report.  The  ALFGEE  coordinate  system,  formed  by  three  Cartesian  coordinates: 
x,y,  and  z,  is  right  handed.  The  x-coordinate,  commonly  called  the  range,  is 

usually  the  distance  measured  from  the  detector  (or  field  point)  where  the 
irradiance  is  being  calculated  to  the  plane  containing  the  nozzle  centers  of 
symmetry.  The  y-coordinate,  commonly  called  the  cross  range,  is  usually  the 
distance  measured  from  the  field  point  to  the  plane  passing  through  the  center 
of  an  even  numbered  array  of  nozzles  that  is  perpendicular  to  the  plane 
containing  the  nozzle  centers  of  symmetry.  Historically,  the  nozzles  have 

been  placed  along  the  ground  in  a  straight  line.  The  z-coordinate,  height,  is 
normally  measured  to  and  above  the  ground  plane. 

Occasionally,  the  origin  of  the  coordinate  sytem  will  be  changed, 

typically  for  an  experimenter's  convenience.  This  is  readily  accomplished, 
since  the  implementation  of  the  algorithms  has  always  included  the  flame 
coordinates. 

The  normal  to  the  detecting  surface  is  shown  in  Figure  31  as  a  unit 

vector,  n.  This  normal  is  described  by  two  angles,  one  is  similar  to  and  is 

called  azimuth  and  the  other  is  its  elevation  above  the  ground  plane.  The 

azimuth  is  specified  by  a  rotation  about  the  height  axis,  where  an  azimuth  of 
zero  is  a  normal  pointed  towards  the  negative  x  direction.  The  elevation  is 
specified  by  a  rotation  about  a  line  perpendicular  to  the  normal  and  in  an  xy 

plane,  ie.  a  plane  of  constant  z. 
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TREATMENT  OF  GROUND  REFLECTIONS  (ALBEDO). 


When  the  flames  are  close  to  other  surfaces,  some  energy  will  be 
scattered  by  these  other  surfaces  and  may  illuminate  the  collecting  surface 
under  consideration.  Usually,  the  TRS  nozzles  are  fired  quite  close  to  the 
ground. 


Calculation  of  the  scattered  radiation  would  be  more  difficult,  if 

done  correctly,  than  the  calculation  of  the  direct  radiation.  This  is  because 

the  ground  is  a  diffuse  reflector.  Radiation  incident  upon  the  ground  is 

scattered  into  the  hemisphere  of  air  above.  This  scattered  radiation  reaches 
the  collecting  surface  if  the  scattering  point  is  visible  to  the  collector. 

Therefore,  the  incident  flux  on  the  ground  should  be  calculated  at  many  points 
on  the  ground,  and  then  an  integration  over  this  large  source  performed  for 
each  collector  location. 

Again  to  reduce  the  amount  of  computation,  certain  approximations 
were  made.  Rather  than  perform  the  integration,  we  decided  to  treat  the 
ground  as  a  specular  reflector  with  a  reflection  coefficient  equal  to  the 
average  albedo.  This  albedo  is  also  a  parameter  that  the  analyst  can 
adjust  so  calibration  will  again  mitigate  some  of  the  errors  associated 
with  this  approach.  However,  this  model  will  not  calculate  any  contribution 
from  the  ground  for  any  collecting  surface  facing  away  from  the  flames.  In 
that  case,  the  model  will  clearly  underpredicf  the  received  flux. 

When  the  decision  was  made  to  proceed  in  this  fashion,  it  was  hoped 
that  some  experimental  data  would  be  available  soon,  to  allow  an  assessment  of 
the  adequacy  of  this  technique.  Unfortunately,  none  of  the  recommended  albedo 
experiments  had  been  made  by  the  time  this  effort  was  completed.  Some 
attempts  at  a  measurement  have  since  been  made  (early  1983),  but  these  are 
inconclusive.  Meaningful  measurements  need  to  be  made.  Nevertheless,  this 
approach  seems  to  admit  to  reasonable  modeling  of  the  data.  At  some  point  when 
more  reliance  is  placed  on  the  modeling,  it  seems  reasonable  to  resolve  the 
remaining  issues. 
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Finally,  this  approach  used  does  not  provide  any  angular  dependence 
of  the  scattering  by  the  ground.  If  the  reflectivity  of  the  ground  changes  as 
the  angle  of  incidence  (off  the  normal)  changes  as  some  authors  have 
suggested,  the  albedo  should  depend  on  ground  range  from  flame  to  collector. 


4-4  STRENGTH  OF  A  NOZZLE. 

The  model  also  includes  another  adjustable  parameter  that  allows  the 
analyst  to  modify  the  power  of  each  flame  and  hence  scale  the  output.  This 
parameter  has  always  been  the  last  one  adjusted  when  calibrating  the  model. 
Once  a  comparison  has  been  made  between  a  subset  of  the  data  base  and  the 
corresponding  calculations  for  one  set  of  parameters,  the  sum  of  the  residuals 
can  always  be  forced  to  be  zero  by  multiplying  the  strength  used  in  the 
calculations  by  the  ratio  of  the  sum  of  the  measured  values  to  the  difference 
of  the  sum  of  the  measured  values  and  the  sum  of  the  residuals. 


4-5  POWER  PER  UNIT  HEIGHT. 

One  of  the  most  important  parts  of  the  model  is  the  distribution  of 
the  radiated  power  per  unit  height  along  the  axis  of  the  flame.  Yet  the 
modeling  to  date  has  relied  on  essentially  no  data  to  establish  this  basic 

relationship.  An  early  cursory  review  of  some  photographic  data  indicated 
that  there  existed  areas  within  the  flame  that  radiated  more  energy  per  unit 
cross  sectional  area  than  other  areas.  Qualitatively,  except  for  a  small 
region  very  close  to  the  nozzles,  the  brightness  in  the  visible  spectrum 
seemed  to  be  a  maximum  near  the  nozzles,  which  decreased  with  height  and  cross 
sectional  width.  Whether  the  flame  was  radiating  as  a  volume  or  as  a  surface 
could  not  be  easily  determined.  However,  if  the  flames  were  indeed  opaque  and 

were  radiating  from  a  relatively  small  region  at  the  surface  of  the  flame,  the 

radiation  clearly  is  not  Lambertian. 


In  spite  of  this  lack  of  data  and  in  spite  of  our  uncertainty  as  to  the 
details  of  how  these  flames  radiate,  we  chose  to  adopt  a  simple  model  based  on 


some  qualitative  interpretations  of  some  photographic  and  video  records.  The 
measure  of  success  in  the  early  mode  ling  efforts  was  the  utility 
demonstrated  by  the  model. 

The  model  of  the  power  per  unit  height  chosen  for  the  cylindrical 
representation  included  two  adjustable  coefficients.  These  coefficients 
provided  a  limited  ability  to  adjust  the  distribution  of  energy  along  the 
flame  axis  to  best  fit  the  available  data.  One  coefficient  was  normalized  by 
the  height  of  the  radiating  flame  and  separated  the  model  into  two  regions.  In 
the  first  region  the  normalized  power  per  unit  length  was  assummed  to  be  a 
constant  value,  unity.  In  the  second  region  the  functional  dependence  chosen 
was  a  linear  decrease  with  height.  The  second  coefficient  was  the  net 
difference  from  unity  of  the  power  per  unit  length  at  a  normalized  height  of 
one.  Thus,  if  the  first  coefficient  equals  one  (or  the  second  equals  zero), 
the  functional  form  would  be  a  constant.  If  the  first  equals  zero,  and  the 
second  is  set  to  one,  the  functional  form  would  be  linear  with  the  maximum  at 
the  nozzle,  decreasing  to  zero  at  the  normalized  flame  height  of  one. 

During  the  development  of  the  inverted  cone  model,  it  was  decided, 
after  much  consideration  of  the  experience  gained  by  calibrating  to  the 
various  data,  that  a  more  realistic  model  should  include  the  width  of  the 
flame  at  each  source  point.  For  the  cylindrical  model  this  of  course  cannot 
make  any  difference.  The  model  that  evolved  was  similar  to  the  original  one 
except  that  the  constant  region  was  removed  and  an  additional  linear  region 
was  added  above  the  radiating  flame.  This  was  accomplished  without  increasing 
the  number  of  coefficients  by  requiring  the  power  per  unit  length  to  decrease 
to  zero  at  the  top  of  the  second  region.  An  alternate  approach  would  have 
been  to  add  a  third  region,  but  this  would  have  required  one  more  adjustable 
parameter  in  the  model. 

4-6  FLAME-FLAME  OBSCURATION. 

The  very  early  models  did  not  reduce  the  received  flux  when  one  flame 
obscurred  the  view  of  another.  If  a  flame  was  completely  shadowed  by  another 


flame,  its  contribution  to  the  flux  was  calculated  as  if  the  closer  flame  were 
transparent  to  the  radiation  from  it.  Initially,  this  approximation  worked 
because  in  use  the  flames  only  infrequently  obscurred  one  another,  or  the 
obscurred  surface  area  was  small.  Eventually,  however,  as  more  nozzles  were 
added  to  the  flame  arrays  and  more  ground  area  was  planned  for  experiments,  it 
became  clear  that  consideration  of  obscuration  effects  was  essential. 

In  the  cylindrical  flame  model,  the  geometry  is  straightforward.  The 
original  algorithms  for  computing  this  geometry  work  well  so  long  as  the 
flames  do  not  touch.  They  lead  to  unphysical  asymmetries  in  the  results  if  the 
flames  touch.  This  early  approach  included  the  assumption  that  the  flames  did 
not  touch,  to  avoid  some  computations.  For  most  of  the  early  cases  this  was 
an  adequate  approximation. 

To  compute  the  radiation  from  any  flame  when  obscuration  effects  are 
included,  the  effect  of  all  the  other  flames  on  it  must  be  considered.  Every 
combination  of  flame  pairs  must  be  evaluated.  In  the  ensuing  discussion,  one 
flame  of  the  pair,  at  which  radiation  is  being  summed,  is  the  radiating  flame 
or  radiator.  The  other  possibly  obscuring  flame  is  called  the  obscurer.  When 
the  flames  do  not  touch,  only  an  obscurer  closer  than  the  radiator  need  be 
included.  Those  farther  out  cannot  obscure  even  if  they  are  different  sizes. 

Rather  than  computing  the  geometry  for  each  combination  of  flame 
pairs,  we  can  avoid  half  of  the  computations  by  sorting  the  flames.  Then  when 
the  range  to  an  obscurer  exceeds  that  for  a  radiator,  in  order  of  range  from 
the  fixed  point  the  rest  of  the  obscurers  for  this  radiator  can  be  ignored. 

When  the  development  of  the  inverted  cone  model  was  started,  the  need 
for  including  intersecting  fiames  was  identified.  Treating  obscuration  is 
more  difficult,  especially  since  the  radius  of  the  flames  depends  upon  height 
in  this  model.  Further,  some  symmetry  is  lost  owing  to  the  ground  albedo. 
These  lead  to  a  substantial  increase  in  computation.  To  reduce  the 
computation,  each  inverted  cone  is  approximated  as  a  set  of  stacked  disks. 
This  simplification  does  not  introduce  appreciable  error. 


It  has  been  assumed  throughout  the  modeling  effort  that  the  models 

would  be  used  for  ensembles  of  identical  flames.  Each  flame  must  be  the  same 
size  and  all  the  nozzles  must  be  at  the  same  height  for  the  obscuration  models 
to  work  correctly.  When  there  exists  a  need  to  use  these  models  for  other 
arrangements  of  TRS  units,  or  mixtures  of  different  types  of  units,  additional 
modeling  effort  will  be  required. 

4-7  AN  ERROR  IN  EARLY  FLAME-FLAME  OBSCURATION  MODEL. 

The  implementation  of  the  ALFGEE  models,  used  for  many  of  the 
calculations  reported  herein,  contains  a  deficiency  in  the  obscuration  model. 

When  flames  touch,  the  obscuration  algorithm  does  r.ot  accurately  calculate  the 
effects.  Table  14  shows  the  results  that  are  generated  with  this 

implementation.  Note  that  for  some  results,  the  expected  symmetry  about  y=0 

is  not  produced.  This  effect  is  seen  in  the  calculations  for  four  nozzles. 
These  calculations  were  performed  for  a  particular  MILL  RACE  experiment. 

This  particular  defect  existed  in  all  obscuration  models  used 
throughout  this  effort.  The  results  presented  in  this  report,  however,  have 
not  been  affected  thereby  with  the  exception  of  the  results  in  Table  14. 

A  listing  of  the  implementation  containing  the  error  is  included  as 
Appendix  E.  This  particular  implementation  was  written  by  the  author  on  an 
Apple  II  Plus  microcomputer.  It  was  written  in  Apple  Pascal  (version  l.l)  for 
use  on  that  machine.  This  listing  is  included  to  complete  the  presentation  of 
the  algorithms  developed  under  this  effort.  These  algorithms  were  implemented 
and  tested  for  actual  cases  of  interest. 


Table  14.  Comparison  of  different  spacings  in  front  of  deckhouse 
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SECTION  5 


TRS  -  AIR  BLAST  INTERFERENCE 

This  section  presents  the  limited  work  that  was  done  to  address  some 
of  the  TRS-air  blast  interference  issues.  This  work  consisted  of  an  initial 
look  at  how  far  out  from  the  TRS  units  blowoff  effects  might  be  seen  on  event 
MILL  RACE.  These  calculations  were  performed  to  estimate  how  much  ground  area 
can  respond  to  the  incident  thermal  radiation.  To  properly  assess  the  effect 
on  the  air  blast  will  require  large  and  long  running  two  and  three  dimensional 
calculations. 

In  addition,  some  simple  one  dimensional  HULL-type  calculations  were 
performed  to  determine  the  magnitude  of  the  effect  a  hot  flame  has  on 
overpressure  when  air  blast  propagates  through  it.  These  calculations  showed 
a  small  effect,  which  is  contrary  to  the  MISERS  BLUFF  result.  This  difference 
is  attributed  to  the  sizable  effect  expected  from  particulates  in  the  flame. 
Unfortunately,  computer  related  problems  precluded  us  from  performing  further 
calculations. 


5-1  BLOWOFF. 

When  thermal  radiation  is  rapidly  deposited  into  the  ground,  it  can 
lead  to  an  effect  known  as  blowoff.  The  deposited  energy  does  not  have 
sufficient  time  to  conduct  away  from  the  irradiated  surface.  For  blowoff  to 
occur,  the  temperature  in  the  ground  must  quickly  increase  to  the  boiling 
point  of  water  and  some  steam  can  be  generated  from  any  free  water. 
Additional  energy  can  raise  the  temperature  further  to  where  water  and  carbon 
dioxide  bound  to  crystalline  material  can  be  explosively  released.  This  rapid 
release  of  gases  can  drag  up  particulates,  which  are  subsequently  heated,  and 
that  appear  as  "blowoff"  material.  The  importance  of  this  material  is  in  how 


it  affects  air  blast  propagating  through  it.  The  nuclear  air  blast  precursor 
is  a  good  example. 


The  whole  blowoff  process  is  extremely  complicated  and  not  well 

understood  even  after  considerable  research.  It  is  known  that  the  amount  of 
energy  that  it  takes  per  unit  area  of  soil  to  cause  blowoff  depends  on  many 

factors  such  as  soil  type,  flux,  and  incident  angle  to  name  a  few. 

Nonetheless,  for  NTS  type  soils,  similar  to  WSMR  soils,  and  for 
nuclear  weapon  flux  levels  (for  tens  of  Kiloton  devices)  an  approximate  value 
of  six  (6)  calories  per  square  centimeter  has  been  used  (Reference  5)  and 

seems  to  be  an  appropriate  value  for  our  purposes  for  estimating  incipient 

blowoff.  This  criterion  is  called  the  blowoff  threshold. 


5-2  POTENTIAL  BLOWOFF  ON  MILL  RACE. 

Figures  32  through  37  present  contours  of  equal  peak  flux  incident  on 
the  ground  calculated  posttest  for  the  MILL  RACE  event.  The  expected  blowoff 
area  is  determined  by  this  peak  flux  multiplied  by  the  burn  time  to  give 
fluence.  Where  the  fluence  exceeds  the  blowoff  threshold,  blowoff  can  occur. 
Each  figure  is  drawn  to  the  same  scale  to  allow  comparison.  The  dimensions 
are  in  meters  and  the  flux  is  in  calories  per  square  centimeter  per  second. 

As  an  example,  Figure  32,  shows  the  flux  perpendicular  to  the  ground 
from  the  BRL  MILL  RACE  four  nozzle  array.  The  let  er  N  shows  the  location  of 
each  nozzle  along  the  y-axis.  If  the  burn  rime  is  one  second  and  the  blowoff 
threshold  is  6  calories  per  square  centimeter,  then  contour  5  encloses  the 
blowoff  area.  Contour  5  crosses  the  x  axis  at  about  3.7  meters  and  crosses 
the  y  axis  at  about  3.8  meters.  The  contours  look  almost  circular  about  the 
origin. 


The  results  for  Figure  32  were  based  on  the  cylindrical  model 
calibrated  to  the  BRL  source.  The  data  were  not  generated  for  negative  y 
values  because  of  the  symmetry  on  either  side  of  the  x  axis. 
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ALPGEE  Cylindrical:  BRL  Model  B:  Flux  perpendicular  to  ground 
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ALFGEE  Cylindrical:  UK1  Model  S3:  Flux  perpendicular  to  ground 
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Figure  53.  Flux  on  ground  contours:  UK- 1  model  .35. 


ALFGEE  Cylindrical:  UK2  Model  C:  Flux  perpendicular 
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Figure  34.  Flux  on  ground  contours:  Uk2  model  C. 
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ALFGEE  Cylindrical:  Navy  Model  C:  Flux  perpendicular  to  ground 
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Figure  55.  Flux  on  /round  contours:  N'avv  model  C. 


ALFGEE  Cylindrical:  Navy  Model  D:  Flux  perpendicular  to  ground 
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ALFGEE  Inverted  Cone:  Navy:  Flux  perpendicular  to  ground 
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Figure  57.  Flux  on  ground  contours:  Navy  inverted  cone. 


Figure  33,  on  the  other  hand,  shows  that  symmetry  did  not  exist 

during  MILL  RACE.  The  UK-!  nozzles  3  and  4  did  not  ignite,  and  nozzle  7 

burned  unreliably.  Adjusting  the  output  for  these  three  nozzles  removes  the 

symmetry,  thus  both  positive  and  negative  y  values  were  computed  and  plotted. 

Knowing  the  burn  time  for  an  array,  the  contour  plots  can  be  used  to 
estimate  the  area  of  blowoff,  given  any  blowoff  threshold  for  burn  times 
similar  to  those  used  in  MILL  RACE.  The  bum  times  for  each  nozzle  array  are: 
Navy  2.4  seconds;  BRL  l.l  seconds;  UK-2  1.0  seconds;  and  UK-1  2.0  seconds 
(Reference  3).  By  way  of  example,  desert  soil  at  NTS  can  start  to  blowoff  at 
about  6  calories  per  square  centimeter.  The  burn  time  for  the  eight  nozzle 
UK-1  MILL  RACE  array  was  2  seconds.  The  above  blowoff  criterion  then 
corresponds  to  a  flux  of  3  calories  per  square  centimeter  per  second,  which  is 
contour  3.  The  extent  along  the  x  axis  is  about  6.6  meters.  The  burn  time, 
on  the  other  hand,  for  the  four  nozzle  UK-2  MILL  RACE  array  was  only  one 
second.  Therefore,  contour  5  represents  the  blowoff  threshold.  It  extends 
out  to  about  4.7  meters. 

The  area  that  exceeds  the  blowoff  threshold  for  UK- 1  greatly  exceeds 
that  for  UK-2  because  the  burn  time  is  twice  as  long  and  there  are  one  to  two 
more  nozzles  burning. 

5-3  ONE  DIMENSIONAL  FLAME  CALCULATIONS. 

Table  15  summarizes  the  cases  run  to  investigate  the  influence  that  a 
slab  of  hot  air,  sandwiched  between  cold  ambient  air,  would  have  on  air  blast 
propagating  through  this  layer. 

Three  incident  overpressures  were  considered:  3,  7,  and  10  psi.  The 
case  number  (first  column  in  the  table)  is  the  overpressure  multiplied  by  100. 
Two  flame  widths  were  briefly  looked  at.  They  were  one  and  two  meters  and  are 
indicated  by  the  second  column.  The  third  column  shows  the  ratio  of  the  flame 
internal  energy  to  ambient.  The  temperature  ratio  is  not  shown  and  depends  on 
the  equation  of  state.  The  corresponding  flame  temperature  is  approximately 
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between  2500  and  3000  degrees  Kelvin,  although  the  results  are  not  too 

sensitive  to  the  temperature  in  this  range. 

The  results  are  shown  for  various  observer  locations  beyond  the 

flame.  Two  blocks  of  data  are  shown  for  each  overpressure.  The  first  (top) 

block  is  the  overpressure  at  the  corresponding  observer  location.  The  second 

block  is  the  percentage  differences.  The  first  line  in  each  block  is  the  free 
air  case,  ie.  no  flame  is  present.  This  is  the  baseline  case  against  which 
every  other  case  is  compared.  The  distance  from  the  center  of  the  flame  to 

observer  110  is  !0. 

The  largest  change  observed  was  for  10  psi  and  this  change  was  only 
three  percent  for  distances  close  to  the  flame.  Thus  the  induced  effect  due 

to  the  presence  of  just  a  hot  air  layer  is  small  for  the  flame  widths 
considered  for  overpressures  less  than  10  psi. 

It  is  our  hypothesis  that  the  TRS  flame  must  have  some  matter  inside 
them  that  does  not  contribute  significantly  to  the  pressure.  Otherwise,  it 

would  be  expected  that  the  bright  hot  flame  would  produce  associated  air 
blast,  as  was  seen  in  early  TRS  bag  designs.  Whether  solid  particulates  are 
formed  or  liquid  droplets  (molten  aluminum  oxide,  for  example)  remains  to  be 

determined.  In  either  case,  it  is  expected  that  once  sufficient  particulates 
or  droplets  have  been  formed  with  a  resultant  material  density  of  the  order  of 
the  air  density,  then  the  flow  will  be  affected  significantly  by  these  small 
pieces  of  matter  through,  at  least,  aerodynamic  drag.  Therefore,  to  establish 
the  non-interference  of  TRS  flames,  additional  one  dimensional  calculations 
with  interactive  particles,  that  can  numerically  simulate  either  solid 
particulates  or  liquid  droplets,  will  need  to  be  performed. 


5-4  ADDITIONAL  AIR  BLAST  INTERFERENCE  WORK. 

It  was  shown  in  the  previous  section  that  the  induced  effect  on  air 
blast  due  to  the  presence  of  just  a  hot  air  layer  is  small  for  the  flame 
widths  considered  for  overpressures  less  than  10  psi.  In  order  to  establish 
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the  non-interference  of  TRS  flames,  additional  one  dimensional  calculations 
with  interactive  particles  will  need  to  be  performed. 

Another  effect  that  needs  to  be  addressed  is  the  induced  three 
dimensional  flow  as  a  result  of  the  three  dimensional  character  of  the  TRS 
flame-air  blast  interaction.  Such  an  investigation  was  outside  of  the  scope  of 
this  effort. 

The  current  practice  in  the  field  is  to  mitigate  an/  air  blast 
interference  problems  by  firing  the  flames  early  enough  to  allow  the  complete 
burn  of  the  flame  as  well  as  its  rise  out  of  the  way  before  air  blast  arrival. 
Although  this  is  felt  to  be  a  reasonable  approach  for  reducing  the  direct 
influence  of  the  flame  on  the  local  peak  overpressure,  the  effect  of  the 
combustion  products  on  air  blast  further  down  range  and  on  the  rest  of  the 
local  air  blast  waveform  has  not  been  established.  The  most  significant 
problem  introduced  by  this  approach  to  air  blast  interference  mitigation  may 
well  be  the  possibly  artificial  delay  introduced  in  the  arrival  of  the  air 
blast  relative  to  thermal  exposure.  Considerable  analysis  may  be  required  for 
each  weapon  system  thusly  exposed  on  a  combined  air  blast  /  thermal  test. 


SECTION  6 


CONCLUSIONS 

This  section  presents  the  conclusions  drawn  in  the  area  of  estimating 
the  radiated  thermal  environments  from  aluminum-liquid  oxygen  (AI-LOX)  thermal 
radiation  simulators  (TRS).  Some  conclusions  were  also  drawn  in  the  area  of 
TRS-air  blast  interference,  however  these  were  significantly  limited. 


6-1  MODEL  DEVELOPMENT. 

Most  importantly,  this  effort  produced  models  to  calculate  flame  TRS 
generated  peak  flux  that  represent  those  same  environments  in  the  field.  The 
inverted  cone  model  is  the  best  developed  so  far  when  flame  obscuration  and  an 
accurate  treatment  of  view  angle  are  included.  This  model  produces  results 
that  are  in  closest  agreement  with  the  data  up  through  MILL  RACE. 

It  is  now  possible  to  use  one  calibration  of  the  inverted  cone  model 
to  represent  all  cases.  Consequently,  the  data  base  may  now  be  of  sufficient 
size  to  allow  certain  statistical  inferences  to  be  made  about  the  inherent 
variability  of  the  flame  TRS  in  use  in  the  field. 

Model  development  was  evolutionary.  First  we  tried  to  provide  a  tool 
for  interpolation  based  on  simple  physical  principles.  Nexr,  refinements  were 
developed  in  the  representation  of  flame  geometry,  and  a  more  accurate 
treatment  of  flame-flame  obscuration.  While  model  development  was  ongoing, 
the  model  was  exercised  for  TRS  development  work,  pre-MILL  RACE  nozzle 
checkout,  experiment  site  layout,  and  experiment  interpolation  and  analysis. 

Some  of  the  early  data  were  affected  to  an  unknown  extent  by 
calorimeter  stand  tip-back.  This  effect  was  more  serious  on  the  late  time 


flux  data  for  much  of  the  early  calorimeter  data.  The  current  practice  is  to 
use  heavier  poles  at  the  TRS  sites;  these  are  much  less  suseptible  to  thermal 
stress  deformation. 

We  conclude  that  a  substantially  more  exhaustive  and  a  more  careful 

effort  to  obtain  accurate  calorimeter  data  needs  to  be  conducted.  Further, 

more  thorough  data  interpretation  and  analysis  is  required.  At  least  four 
different  individuals  have  read  the  "peak"  flux  off  of  oscilliscope  traces 
that  have  been  entered  into  the  DNA  data  base.  Now  that  these  data  are  being 
recorded  digitally,  the  detailed  temporal  waveforms  can  be  made  available  for 
more  meaningful  data  analysis.  There  has  been  no  consistent  approach  applied 
to  evaluation  of  the  available  data.  In  some  cases  data  was  withheld  as  not 
being  representative. 

The  TRS  development  to  date  has  been  directed  toward  improving  the 
flame  generation.  Little  emphasis  has  been  applied  toward  evaluating 

reliability  or  repeatability.  We  strongly  recommend  that  one  design  be  fixed 

for  future  use  in  the  DNA  HE  program,  and  that  this  design  be  calibrated  under 
realistic  field  conditions.  Such  calibration  should  be  done  using  many  more 
calorimeters  than  have  been  used  in  the  past.  While  it  is  recognized  that 
funding  limitations  tend  to  preclude  having  enough  of  what  is  needed,  the 
calibration  of  the  DNA  thermal  simulation  "workhorse"  should  be  viewed  as  the 
foundation  of  the  nuclear  thermal  effects  simulation  program. 

6-2  TRS-AIR  BLAST  INTERFERENCE. 

In  the  area  of  TRS-Air  Blast  interference,  the  results  are  less 
clear.  The  conducted  effort  is  largely  incomplete,  and  additional  work  is 
required  to  resolve  the  issues.  Calculations  indicate  *hat  incipient  blowoff 
of  ground  material  was  probable  on  MILL  RACE.  Other  calculations  have 
indicated  that  incipient  ground  blowoff  is  probable  for  the  normal  one  second 
burn  at  the  DNA  TRS  site  at  KAFB.  Additional  burn  time  would  extend  the 
blowoff  area.  What  looks  like  the  "steam  porch"  phase  of  ground  blowoff  has 
been  seen  in  the  photographs  taken  at  the  TRS  site. 


It  was  also  concluded  that  modeling  a  TRS  flame  as  merely  a  hot  clean 
flame,  i.e.  one  containing  no  particulates,  will  not  accurately  simulate 
results  from  MISERS  BLUFF.  During  MISERS  BLUFF  the  overpressure  was  reduced 
on  one  of  the  targets  placed  immediately  behind  the  TRS  units.  Numerical 
calculations  performed  for  overpressures  up  to  10  psi  and  for  flame  widths  up 
to  two  meters  wide  demonstrated  that  other  effects  will  need  to  be  considered. 
Additional  work  is  recommended  to  quantitatively  understand  the  influence  of 
particulates  on  the  effect  seen  during  MISERS  BLUFF. 

Finally,  much  difficulty  in  accessing  the  DNA-provided  computational 
facilities  at  KAFB  made  it  very  difficult  to  continue  the  investigation  of 
TRS-air  blast  interference. 
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SECTION  7 


RECOMMENDATIONS 

This  section  presents  recommendations  that  resulted  from  performance 
rf  the  work. 

We  recommend  that  additional  work  be  performed  to  quantitatively 
understand  the  influence  of  particulates  in  the  effect  seen  during  MISERS 
BLUFF. 


Our  recommendations  in  the  area  of  TRS  modeling  is  that  much  more 
experimental  and  some  modeling  work  needs  to  be  done.  In  the  area  of  TRS- 
airblast  interference,  our  recommendation  is  to  determine  how  well  the  current 
interference  mitigation  techniques  will  find  acceptance  in  the  nuclear 
survivability  and  vulnerability  community. 

A  substantially  more  exhaustive  and  a  more  careful  effort  to  obtain 
calorimeter  data  needs  to  be  done.  This  recommendation  is  made  with  the 
knowledge  that  the  data  base  essentially  doubled  during  September  1982  when 
DNA  and  SAI  conducted  a  Calibration  Program  at  the  DNA  TRS  site  at  Kirtland 
AF8.  Even  with  those  data,  significant  additional  data  under  more  controlled 
conditions  is  necessary  if  the  community  expects  to  place  great  reliance  on 
modeled  results. 

More  thorough  data  interpretation  and  analysis  is  recommended.  At 
least  four  different  individuals  have  read  the  "peak"  flux  off  of  oscilliscope 
traces  that  have  been  entered  into  the  DNA  data  base.  Now  that  these  data  are 
being  recorded  digitally,  the  detailed  temporal  waveforms  could  be  made 
available  for  more  meaningful  data  analysis.  There  has  been  no  consistent 
approach  applied  to  evaluating  these  data. 


108 


The  TRS  development  to  date  has  been  geared  to  improving  the  flame 
generation.  Little  emphasis  has  been  given  to  evaluating  reliability  or 
repeatability.  We  strongly  recommend  that  a  design  be  frozen  for  future  use 
on  the  DNA  HE  program,  and  that  that  design  be  calibrated  under  realistic 
field  conditions. 

Finally,  we  regretably  observe  that  much  difficulty  in  accessing  the 
DNA-provided  computational  facilities  at  KAFB  made  it  virtually  impossible  to 
continue  the  investigation  of  TRS-air  blast  interference.  It  is  recommended 
that  DNA,  not  a  contractor,  more  closely  monitor  difficulties  that  contractors 
are  having  when  using  the  government  furnished  equipment.  Whether  this  would 
lead  to  any  improvement  in  throughput  is  debatable;  nevertheless,  we  feel  it 
is  important  to  improve  this  communication  of  the  problems  we  have  been 
having,  so  that  DNA  management  might  have  better  information  available  when 
decisions  are  made. 
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APPENDIX  A 


ANALYSIS  OF  CALORIMETER  POLE  BENDING 

This  appendix  presents  results  from  an  analysis  of  the  influence  of 

the  bending  of  the  pipe,  that  holds  the  fixtures  that  in  turn  hold  the 
calorimeters,  on  received  thermal  flux.  The  bending  of  a  pipe,  commonly 

called  a  calorimeter  stand,  is  due  to  thermal  stress  that  results  from 
asymmetrical  thermal  loads.  The  effects  can  be  significant.  For  this  reason, 
the  calorimeter  stands  were  replaced  with  larger  pipe,  which  will  not  be  so 
easily  bent. 

The  analysis  consisted  of  inspecting  the  sensitivity  of  the  flux  to 

the  changes  of  a  pole's  position  as  a  result  of  irradiation.  To  simplify  the 
analysis,  the  case  considered  was  for  uniform  illumination  of  one  side  of  the 

pole,  causing  it  to  tip  back  with  a  shape  identical  to  an  arc  on  a  circle. 
This  shape  change  causes  the  calorimeter's  range  to  increase,  its  height  to 
decrease,  and  its  normal's  elevation  to  increase. 

The  result  is  that  each  of  the  three  above  changes  will  reduce  the 

received  flux  as  the  bending  increases.  The  dominant  effect,  considering  the 

change  in  geometry  relationships,  is  the  effect  on  range.  For  every 
centimeter  the  pole  is  tipped  back  at  the  calorimeter,  the  flux  was  reduced 
0.4%.  During  the  early  testing,  we  estimate  that  this  effect  could  have 
been  of  the  order  of  10  percent  for  later  times.  This  is  another  reason 

for  choosing  to  model  the  early  peak  flux,  rather  than  the  fluence.  Recent 
data  are  not  affected  by  this  phenomenon  due  to  larger  pipe  size.  Table  16 

summarizes  the  type  of  results  generated. 
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THJi  TEMPLATE  IS  USED  TO  CALCULATE  EFFECT  OF  BEAM  BENDING  ON  PEAK  FLUX. 


APPENDIX  B 


TRS  PEAK  FLUX  DATA:  PRIOR  TO  MILL  RACE 

This  appendix  presents  peak  flux  data  taken  on  various  two  and  four 

nozzles  burns  performed  on  the  current  DNA  TRS  nozzles  at  the  KAFB  site  prior 
to  the  testing  of  the  MILL  RACE  nozzles.  This  appendix  documents  where  each 
calorimeter  was  located,  as  well  as  different  interpretations  of  the  "peak" 
flux's  value. 

Three  tables  are  included.  The  first,  Table  17,  documents  events  D- 

15  through  D-45,  where  data  of  some  significance  existed.  Comparisons  were 
made  with  an  early  calibration  of  the  right  cylinder  model  in  this  table. 

During  the  period  of  time  when  these  data  were  being  recorded,  numerous 
adjustments  were  made  to  the  operating  parameters  for  the  nozzles,  and 
therefore,  the  variability  in  the  data  is  considerable. 

The  second  table,  Table  18,  presents  the  data  taken  during  the 
Pentolite  series.  Comparisons  were  also  made  with  an  early  calibration  of  the 
right  cylinder  model  in  this  table.  These  data  are  interesting  because  more 

calorimeters  were  available  on  each  shot. 

Finally,  Table  19,  re-documents  some  of  the  events  in  the  series  D-34 
through  D-45.  However,  here  comparisons  were  made  with  an  early  calibration 
of  the  inverted  cone  model.  On  the  average,  the  inverted  cone  model  improved 
agreement  with  the  data. 
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APPENDIX  C 


A  MODEL  COMPARED  WITH  DATA  BASE 

This  appendix  presents  a  comparison  of  the  available  peak  flux  data 
with  one  calibration  to  the  inverted  cone  model.  Two  tables  are  included; 
both  list  the  same  data,  but  in  two  different  orders.  One  is  sorted  by  event, 
Table  20,  and  is  almost  chronological.  The  other,  Table  21,  is  sorted  by  the 
ratio  of  model  results  to  measured  results.  The  model  used  was  the  inverted 
cone.  These  lists  illustrate  how  the  differences  are  statistically  dispersed. 


Table  20.  Event  sort  of  model  comparison  with  data  base 


X 


R  CALOR  EVENT 


-182.90 

152.40 

1.11 

FX-1 

D15 

365.80 

304.80 

0.54 

FX-2 

Dl5 

182.90 

304.80 

0.69 

FX-3 

Dl5 

182.90 

152.40 

1.09 

FX-4 

D15 

-182.90 

152.40 

1.21 

FX-1 

D17 

365.80 

304.80 

0.82 

FX-2 

Dl7 

182.90 

304.80 

0.76 

FX-3 

Dl7 

182.90 

152.40 

1.13 

FX-4 

D17 

-182.90 

152.40 

1.16 

FX-1 

D19 

365.80 

304.80 

0.75 

FX-2 

D19 

182.90 

304.80 

0.78 

FX-3 

Dl9 

182.90 

152.40 

1.13 

FX-4 

D19 

-182.90 

152.40 

1.16 

FX-1 

D20 

365.80 

304.80 

1.00 

FX-2 

D20 

182.90 

304.80 

0.88 

FX-3 

D20 

182.90 

152.40 

1.26 

FX-4 

D20 

-182.90 

152.40 

1.21 

FX-1 

D21 

365.80 

304.80 

0.80 

FX-2 

D21 

182.90 

304.80 

0.78 

FX-3 

D21 

182.90 

152.40 

1.17 

FX-4 

D21 

121.90 

228.60 

1.11 

FX-1 

D23 

365.80 

304.80 

0.86 

FX-2 

D23 

182.90 

304.80 

0.83 

FX-3 

D23 

182.90 

152.40 

1.22 

FX-4 

D23 

121.90 

228.60 

1.27 

FX-1 

D24 

365.80 

304.80 

0.96 

FX-2 

D24 

182.90 

304.80 

0.97 

FX-3 

D24 

182.90 

152.40 

1.26 

FX-4 

D24 

121.90 

228.60 

1.11 

FX-1 

D25 

365.80 

304.80 

0.91 

FX-2 

D25 

182.90 

304.80 

0.83 

FX-3 

D25 

182.90 

152.40 

1.05 

FX-4 

D25 

121.90 

228.60 

1.11 

FX-1 

D26 

365.80 

304.80 

0.89 

FX-2 

D26 

182.90 

304.80 

0.85 

FX-3 

D26 

182.90 

152.40 

1.05 

FX-4 

D26 

121.90 

228.60 

1.11 

FX-1 

D28 

365.80 

304.80 

0.97 

FX-2 

D28 

182.90 

304.80 

0.93 

FX-3 

D28 

182.90 

152.40 

1.26 

FX-4 

D28 

121.90 

228.60 

1.11 

FX-1 

D29 

365.80 

304.80 

0.86 

FX-2 

D29 

182.90 

304.80 

0.83 

FX-3 

D29 

182.90 

152.40 

1.17 

FX-4 

D29 

121.90 

228.60 

1.11 

FX-1 

D30 

365.80 

304.80 

0.81 

FX-2 

D30 

182.90 

304.80 

0.88 

FX-3 

D30 

182.90 

152.40 

1.17 

FX-4 

D30 

121.90 

228.60 

1.27 

FX-1 

D3 1 

365.80 

304.80 

1.08 

FX-2 

D31 

182.90 

304.80 

1.10 

FX-3 

D31 

182.90 

152.40 

1.00 

FX-4 

D31 

Table  20.  Event  sort  of  model  comparison  with  data  base  (continued) 


121.90 

228.60 

1.30 

FX-1 

D32 

365.80 

304.80 

1.09 

FX-2 

D32 

182.90 

304.80 

1.11 

FX-3 

D32 

182.90 

152.40 

1.22 

FX-4 

D32 

121.90 

228.60 

1.19 

FX-1 

D33 

365.80 

304.80 

1.08 

FX-2 

D33 

182.90 

304.80 

0.97 

FX-3 

D33 

182.90 

152.40 

1.17 

FX-4 

D33 

365.80 

228.60 

1.18 

FX-3 

D35 

365.80 

228.60 

1.07 

FX-4 

D35 

152.40 

228.60 

1.24 

FX-1 

D34 

365.80 

304.50 

1.25 

FX-2 

D34 

182.90 

304.50 

1.05 

FX-3 

D34 

182.90 

152.40 

1.25 

FX-4 

D34 

365.80 

228.60 

1.16 

FX-1 

D39 

365.80 

228.60 

1.18 

FX-2 

D39 

365.80 

228.60 

1.18 

FX-3 

D39 

365.80 

228.60 

1.02 

FX-4 

D39 

182.90 

228.60 

1.00 

FX-1 

D40 

182.90 

228.60 

1.06 

FX-2 

D40 

182.90 

228.60 

0.86 

FX-3 

D40 

182.90 

228.60 

0.74 

FX-4 

D40 

274.30 

228.60 

1.30 

FX-1 

D43 

274.30 

228.60 

1.21 

FX-2 

D43 

274.30 

228.60 

0.91 

FX-3 

D43 

274.30 

228.60 

0.73 

FX-4 

D43 

304.80 

228.60 

1.06 

FX-1 

D44 

182.90 

203.00 

1.19 

FX-3 

D44 

304.80 

228.60 

0.99 

FX-1 

D45 

182.90 

203.00 

1.13 

FX-4 

D45 

200. 00 

200.00 

1.09 

FX-1 

PENT- 2 

200.00 

200.00 

0.95 

FX-2 

PENT- 2 

200.00 

200.00 

1.15 

FX-3 

PENT- 2 

200.00 

200.00 

1.09 

FX-4 

PENT- 2 

350.00 

200.00 

0.95 

FX-5 

PENT- 2 

350.00 

200.00 

0.97 

FX-6 

PENT-2 

350.00 

200.00 

1.04 

FX-7 

PENT-2 

500.00 

200.00 

0.92 

FX-11 

PENT-2 

200.00 

200.00 

0.95 

FX-1 

PENT-3 

200.00 

200.00 

1.08 

FX-2 

PENT-3 

200.00 

200.00 

0.96 

FX-3 

PENT-3 

200.00 

200.00 

0.66 

FX-4 

PENT-3 

350.00 

200.00 

1.06 

FX-5 

PENT-3 

350.00 

200.00 

1.01 

FX-6 

PENT-3 

350.00 

200.00 

0.91 

FX-7 

PENT-3 

500.00 

200.00 

1.04 

FX-8 

PENT-3 

500.00 

200.00 

1.09 

FX-10 

PENT-3 

500.00 

200.00 

1.04 

FX-11 

PENT-3 

200.00 

200.00 

0.77 

FX-1 

PENT-4 

200.00 

200.00 

1.12 

FX-2 

PENT- 4 

200.00 

200.00 

1.22 

FX-3 

PENT-4 

200.00 

200.00 

0.77 

FX-4 

PENT-4 

350.00 

200.00 

0.98 

FX-5 

PENT-4 

350.00 

200.00 

1.22 

FX-6 

PENT-4 

350.00 

200.00 

1.06 

FX-7 

PENT-4 

500.00 

200.00 

1.04 

FX-8 

PENT-4 

500.00 

200.00 

1.38 

FX-10 

PENT-4 

Table  20.  Event  sort  of  model  comparison  with  data  base  (continued) 


500.00 

200.00 

1.14 

FX-11 

PENT-4 

200.00 

200.00 

1.01 

FX-1 

PENT-5 

200.00 

200.00 

1.01 

FX-2 

PENT- 5 

200.00 

200.00 

1.22 

FX-3 

PENT-5 

200.00 

200.00 

1.13 

FX-4 

PENT-5 

350.00 

200.00 

1.21 

FX-5 

PENT-5 

350.00 

200.00 

1.17 

FX-6 

PENT-5 

350.00 

200.00 

1.29 

FX-7 

PENT-5 

500.00 

200.00 

1.04 

FX-8 

PENT-5 

500.00 

200.00 

1.16 

FX-10 

PENT-5 

500.00 

200.00 

1.25 

FX-11 

PENT-5 

300.00 

317.50 

1.22 

FX-1 

UKl-PRE 

300.00 

317.50 

1.26 

FX-2 

UK1-PRE 

300.00 

317.50 

1.25 

FX-3 

UKl-PRE 

269.00 

317.50 

1.30 

CAL -4 

UKl-PRE 

183.00 

317.50 

1.14 

CAL-6 

UKl-PRE 

300.00 

317.50 

1.03 

FX-4 

UKl-PRE 

589.00 

317.50 

1.24 

CAL-1 

UKl-PRE 

511.00 

317.50 

1.20 

CAL-2 

UKl-PRE 

511.00 

317.50 

1.28 

CAL-3 

UKl-PRE 

589.00 

317.50 

1.29 

CAL-5 

UKl-PRE 

0.00 

317.50 

0.57 

CAL-1 

NAV-PRE 

244.00 

317.50 

1.05 

CAL-2 

NAV-PRE 

589.00 

317.50 

1.08 

CAL-3 

NAV-PRE 

244.00 

317.50 

1.15 

CAL-4 

NAV-PRE 

244.00 

317.50 

1.08 

CAL-5 

NAV-PRE 

244.00 

317.50 

0.81 

CAL-6 

NAV-PRE 

244.00 

317.50 

1.07 

FX-1 

NAV-PRE 

589.00 

317.50 

1.07 

FX-2 

NAV-PRE 

589.00 

317.50 

1.12 

FX-3 

NAV-PRE 

244.00 

317.50 

1.00 

FX-4 

NAV-PRE 

274.00 

335.00 

0.95 

CAL-1 

UK2-PRE 

274.00 

122.00 

0.99 

CAL-2 

UK2-PRE 

274.00 

335.00 

1.01 

CAL-3 

UK2-PRE 

274.00 

122.00 

1.05 

CAL-4 

UK2-PRE 

274.00 

335.00 

1.01 

CAL-5 

UK2-PRE 

274.00 

122.00 

0.96 

CAL-6 

UK2-PRE 

274.00 

122.00 

0.89 

CAL-7 

UK2-PRE 

274.00 

229.00 

0.88 

FX-1 

UK2-PRE 

274.00 

229.00 

0.93 

FX-2 

UK2-PRE 

274.00 

229.00 

0.95 

FX-3 

UK2-PRE 

274.00 

229.00 

1.03 

FX-4 

UK2-PRE 

183.00 

172.00 

1.22 

CAL-1 

BRL-PRE 

183.00 

265.00 

1.09 

CAL-2 

BRL-PRE 

183.00 

357.00 

1.08 

CAL-3 

BRL-PRE 

183.00 

357.00 

1.20 

CAL-4 

BRL-PRE 

183.00 

357.00 

1.13 

CAL-5 

BRL-PRE 

183.00 

265.00 

1.04 

CAL-6 

BRL-PRE 

183.00 

172.00 

1.11 

CAL-7 

BRL-PRE 

183.00 

265.00 

1.06 

FX-1 

BRL-PRE 

183.00 

200.00 

1.22 

FX-2 

BRL-PRE 

183.00 

172.00 

1.16 

FX-3 

BRL-PRE 

183.00 

200.00 

1.17 

FX-4 

BRL-PRE 

180.30 

116 . 20 

0.82 

F-l 

UK-1  MR 

271.80 

166.40 

0.94 

F-2 

UK-1  MR 

261.60 

139.10 

1.03 

F-3 

UK-1  MR 

596.90 

229.90 

0.81 

F-4 

UK-1  MR 

Table  20.  Event  sort  of  model  comparison  with  data  base  (continued) 


591.80 

230.50 

0.81 

F-5 

UK-1 

MR 

467.40 

168.90 

1.15 

F-6 

UK-1 

MR 

576.60 

229.90 

1.36 

F-7 

UK-1 

MR 

576.60 

229.20 

1.42 

F-8 

UK-1 

MR 

268.00 

179.10 

0.83 

F-l 

UK-2 

MR 

269.20 

179.10 

0.98 

F-2 

UK-2 

MR 

269.90 

189.90 

1.00 

F-3 

UK-2 

MR 

259.10 

179.10 

0.94 

F-4 

UK-2 

MR 

-190.50 

200.70 

1.67 

F-l 

BRL 

MR 

-190.50 

200.70 

1.07 

F-2 

BRL 

MR 

-190.50 

200.70 

1.00 

F-3 

BRL 

MR 

276.90 

92.70 

0.87 

F-l 

NAVY 

MR 

576.60 

207.00 

1.11 

F-2 

NAVY 

MR 

576.60 

171.40 

1.04 

F-3 

NAVY 

MR 

271.80 

92.70 

0.65 

F-4 

NAVY 

MR 

Table  21.  Ratio 


X  2 

365.80  304.80 

0.00  317.50 

271.80  92.70 

200.00  200.00 

182.90  304.80 

274.30  228.60 

182.90  228.60 

365.80  304.80 

182.90  304.80 

200.00  200.00 
200.00  200.00 

182.90  304.80 

182.90  304.80 

365.80  304.80 

365.80  304.80 

244.00  317.50 

591.80  230.50 

596.90  229.90 

365.80  304.80 

180.30  116.20 

182.90  304.80 

182.90  304.80 

182.90  304.80 

268.00  179.10 

182.90  304.80 

365.80  304.80 

365.80  304.80 

182.90  228.60 

276.90  92.70 

182.90  304.80 

182.90  304.80 

274.00  229.00 

365.80  304.80 

274.00  122.00 

365.80  304.80 

274.30  228.60 

350.00  200.00 

500.00  200.00 

182.90  304.80 

274.00  229.00 

271.80  166.40 

259.10  179.10 

200.00  200.00 

350.00  200.00 

200.00  200.00 

274.00  335.00 

274.00  229.00 

365.80  304.80 

200.00  200.00 

274.00  122.00 

182.90  304.80 


sort  of  model 


R  CALOR 


0.54 

FX-2 

0.57 

CAL-1 

0.65 

F-4 

0.66 

FX-4 

0.69 

FX-3 

0.73 

FX-4 

0.74 

FX-4 

0.75 

FX-2 

0.76 

FX-3 

0.77 

FX-1 

0.77 

FX-4 

0.78 

FX-3 

0.78 

FX-3 

0.80 

FX-2 

0.81 

FX-2 

0.81 

CAL-6 

0.81 

F-5 

0.81 

F-4 

0.82 

FX-2 

0.82 

F-l 

0.83 

FX-3 

0.83 

FX-3 

0.83 

FX-3 

0.83 

F-l 

0.85 

FX-3 

0.86 

FX-2 

0.86 

FX-2 

0.86 

FX-3 

0.87 

F-l 

0.88 

FX-3 

0.88 

FX-3 

0.88 

FX-1 

0.89 

FX-2 

0.89 

CAL-7 

0.91 

FX-2 

0.91 

FX-3 

0.91 

FX-7 

0.92 

FX-11 

0.93 

FX-3 

0.93 

FX-2 

0.94 

F-2 

0.94 

F-4 

0.95 

FX-2 

0.95 

FX-5 

0.95 

FX-1 

0.95 

CAL-1 

0.95 

FX-3 

0.96 

FX-2 

0.96 

FX-3 

0.96 

CAL-6 

0.97 

FX-3 

ison  with  data  base 


EVENT 

D15 
NAV-PRE 
NAVY  MR 
PENT- 3 
D15 
D43 
D40 
D19 
D17 
PENT-4 
PENT-4 
D21 
D19 
D21 
D30 
NAV-PRE 
UK-1  MR 
UK-1  MR 
D17 
UK-1  MR 
D25 
D23 
D29 
UK-2  MR 
D26 
D23 
D29 
D40 
NAVY  MR 
D20 
D30 
UK2-PRE 
D26 
UK2-PRE 
D25 
D4  3 
PENT- 3 
PENT- 2 
D28 
UK2-PRE 
UK-1  MR 
UK-2  MR 
PENT- 2 
PENT- 2 
PENT- 3 
UK2-PRE 
UK2-PRE 
D24 
PENT-3 
UK2-PRE 
D24 


Table  21.  Ratio  sort  of  model  comparison  with  data  base  (continued) 


365.80 

304.80 

0. 

182.90 

304.80 

0. 

350.00 

200.00 

0. 

350.00 

200.00 

0  . 

269.20 

179.10 

0. 

304.80 

228.60 

0. 

274.00 

122.00 

0. 

365.80 

304.80 

1. 

182.90 

152.40 

1. 

182.90 

228.60 

1. 

244.00 

317.50 

1. 

269.90 

189.90 

1. 

-190.50 

200.70 

1. 

350.00 

200.00 

1. 

200.00 

200.00 

1. 

200.00 

200.00 

1. 

274.00 

335.00 

1. 

274.00 

335.00 

1. 

365.80 

228.60 

1. 

300.00 

317.50 

1. 

274.00 

229.00 

1. 

261.60 

139.10 

1. 

350.00 

200.00 

1. 

500.00 

200.00 

1. 

500.00 

200.00 

1 

500.00 

200.00 

1 

500.00 

200.00 

1 

183.00 

265.00 

1 

576.60 

171.40 

1 

182.90 

152.40 

1 

182.90 

152.40 

1 

182.90 

304.50 

1 

244.00 

317.50 

1 

274.00 

122.00 

1 

182.90 

228.60 

1 

304.80 

228.60 

1 

350.00 

200.00 

1 

350.00 

200.00 

1 

183.00 

265.00 

1 

365.80 

228.60 

1 

244.00 

317.50 

1 

589.00 

317.50 

1 

-190.50 

200.70 

1 

365.80 

304.80 

1 

365.80 

304.80 

1 

200.00 

200.00 

1 

589.00 

317.50 

1 

244.00 

317.50 

1 

183.00 

357.00 

1 

182.90 

152.40 

1 

365.80 

304.80 

1 

200.00 

200.00 

1 

200.00 

200.00 

1 

500.00 

200.00 

1 

183.00 

265.00 

1 

182.90 

304.80 

1 

-182.90 

152.40 

1 

97 

FX-2 

D  ZV 

97 

FX-3 

D33 

97 

FX-6 

PENT- 2 

98 

FX-5 

PENT- 4 

98 

F-2 

UK-2  MR 

99 

FX-1 

D45 

99 

CAL-2 

UK2-PRE 

00 

FX-2 

D20 

00 

FX-4 

D31 

00 

FX-1 

D40 

00 

FX-4 

NAV-PRE 

00 

F-3 

UK-2  MR 

00 

F-3 

BRL  MR 

01 

FX-6 

PENT-3 

01 

FX-1 

PENT- 5 

01 

FX-2 

PENT- 5 

01 

CAL-3 

UK2-PRE 

01 

CAL-5 

UK2-PRE 

02 

FX-4 

D39 

03 

FX-4 

UK1-PRE 

03 

FX-4 

UK2-PRE 

03 

F-3 

UK-1  MR 

04 

FX-7 

PENT- 2 

04 

FX-8 

PENT-3 

04 

FX-11 

PENT-3 

04 

FX-8 

PENT-4 

04 

FX-8 

PENT-5 

04 

CAL-6 

BRL-PRE 

.04 

F-3 

NAVY  MR 

.05 

FX-4 

D25 

.05 

FX-4 

D26 

.05 

FX-3 

D34 

.05 

CAL-2 

NAV-PRE 

.05 

CAL-4 

UK2-PRE 

.06 

FX-2 

D40 

.06 

FX-1 

D44 

.06 

FX-5 

PENT-3 

.06 

FX-7 

PENT-4 

.06 

FX-1 

BRL-PRE 

.07 

FX-4 

D35 

.07 

FX-1 

NAV-PRE 

.07 

FX-2 

NAV-PRE 

.07 

F-2 

BRL  MR 

.08 

FX-2 

D31 

.08 

FX-2 

D33 

.08 

FX-2 

PENT-3 

.08 

CAL-3 

NAV-PRE 

.08 

CAL-5 

NAV-PRE 

.08 

CAL-3 

BRL-PRE 

.09 

FX-4 

D15 

.09 

FX-2 

D32 

.09 

FX-1 

PENT-2 

.09 

FX-4 

PENT- 2 

.09 

FX-10 

PENT-3 

.09 

CAL-2 

BRL-PRE 

.10 

FX-3 

D31 

.  11 

FX-1 

Dl  5 

Table  21.  Ratio  sort  of  model  comparison  with  data  base  (continued) 


121.90 

228.60 

1.11 

FX-1 

D23 

121.90 

228.60 

1.11 

FX-1 

D25 

121.90 

228.60 

1.11 

FX-1 

D26 

121.90 

228.60 

1.11 

FX-1 

D28 

121.90 

228.60 

1.11 

FX-1 

D29 

121.90 

228.60 

1.11 

FX-1 

D30 

182.90 

304.80 

1.11 

FX-3 

D32 

183.00 

172.00 

1.11 

CAL-7 

BRL-PRE 

576.60 

207.00 

1.11 

F-2 

NAVY  MR 

200.00 

200.00 

1.12 

FX-2 

PENT-4 

589.00 

317.50 

1.12 

FX-3 

NAV-PRE 

182.90 

152.40 

1.13 

FX-4 

D17 

182.90 

152.40 

1.13 

FX-4 

D19 

182.90 

203.00 

1.13 

FX-4 

D45 

200.00 

200.00 

1.13 

FX-4 

PENT- 5 

183.00 

357.00 

1.13 

CAL-5 

BRL-PRE 

500.00 

200.00 

1.14 

FX-11 

PENT-4 

183.00 

317.50 

1.14 

CAL -6 

UKl-PRE 

200.00 

200.00 

1.15 

FX-3 

PENT- 2 

244.00 

317.50 

1.15 

CAL-4 

NAV-PRE 

467.40 

168.90 

1.15 

F-6 

UK-1  MR 

-182.90 

152.40 

1.16 

FX-1 

D19 

-182.90 

152.40 

1.16 

FX-1 

D20 

365.80 

228.60 

1.16 

FX-1 

D39 

500.00 

200.00 

1.16 

FX-10 

PENT- 5 

183.00 

172.00 

1.16 

FX-3 

BRL-PRE 

182.90 

152.40 

1.17 

FX-4 

D21 

182.90 

152.40 

1.17 

FX-4 

D29 

182.90 

152.40 

1.17 

FX-4 

D30 

182.90 

152.40 

1.17 

FX-4 

D33 

350.00 

200.00 

1.17 

FX-6 

PENT- 5 

183.00 

200.00 

1.17 

FX-4 

BRL-PRE 

365.80 

228.60 

1.18 

FX-3 

D35 

365.80 

228.60 

1.18 

FX-2 

D39 

365.80 

228.60 

1.18 

FX-3 

D39 

121.90 

228.60 

1.19 

FX-1 

D33 

182.90 

203.00 

1.19 

FX-3 

D44 

511.00 

317.50 

1.20 

CAL-2 

UKl-PRE 

183.00 

357.00 

1.20 

CAL-4 

BRL-PRE 

-182.90 

152.40 

1.21 

FX-1 

D21 

-182.90 

152.40 

1.21 

FX-1 

D17 

274.30 

228.60 

1.21 

FX-2 

D43 

350.00 

200.00 

1.21 

FX-5 

PENT- 5 

182.90 

152.40 

1.22 

FX-4 

D23 

182.90 

152.40 

1.22 

FX-4 

D32 

200.00 

200.00 

1.22 

FX-3 

PENT-4 

350.00 

200.00 

1.22 

FX-6 

PENT-4 

200.00 

200.00 

1.22 

FX-3 

PENT-5 

300.00 

317.50 

1.22 

FX-1 

UKl-PRE 

183.00 

172.00 

1.22 

CAL- 1 

BRL-PRE 

183.00 

200.00 

1.22 

FX-2 

BRL-PRE 

152.40 

228.60 

1.24 

FX-1 

D34 

589.00 

317.50 

1.24 

CAL-1 

UKl-PRE 

365.80 

304.50 

1.25 

FX-2 

D34 

182.90 

152.40 

1.25 

FX-4 

D34 

500.00 

200.00 

1.25 

FX-11 

PENT-5 

300.00 

317.50 

1.25 

FX-3 

UKl-PRE 

Table  21.  Ratio  sort  of  model  comparison  with  data  base  (continued) 


182.90 

152.40 

1.26 

FX-4 

182.90 

152.40 

1.26 

FX-4 

182.90 

152.40 

1.26 

FX-4 

300.00 

317.50 

1.26 

FX-2 

121.90 

228.60 

1.27 

FX-1 

121.90 

228.60 

1.27 

FX-1 

511.00 

317.50 

1.28 

CAL-3 

350.00 

200.00 

1.29 

FX-7 

589.00 

317.50 

1.29 

CAL- 5 

121.90 

228.60 

1.30 

FX-1 

274.30 

228.60 

1.30 

FX-1 

269.00 

317.50 

1.30 

CAL-4 

576.60 

229.90 

1.36 

F-7 

500.00 

200. 00 

1.38 

FX-10 

576.60 

229.20 

1.42 

F-8 

190.50 

200.70 

1.67 

F-l 

D20 

D24 

D28 

UK1-PRE 

D24 

D31 

UKi-PRE 

PENT- 5 

.*  *• 

UKI-PRE 

D32 

D43 

UKI-PRE 

UK-1  MR 

PENT-4 

UK-1  MR 

BRL  MR 

APPENDIX  D 


RALE. 

nozzles 


TRS  CALORIMETER  LOCATIONS:  MILL  RACE 

This  appendix  documents  the  locations  of  the  calorimeters  on  MILL 
Four  arrays  of  nozzles  were  fielded.  Three  arrays  consisted  of  four 
each,  and  one  array  consisted  of  eight  nozzles. 
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APPENDIX  E 


ALFGEE  LISTINGS 

This  appendix  presents  listings  of  an  implementation  of  the  ALFGEE 

algorithms  on  an  Apple  II  microcomputer  system.  The  language  used  is  Pascal 
and  the  operating  system  is  Apple  Pascal  Version  l.l.  These  are  derivatives 
of  the  UCSD  II.  I  Pascal  system. 

The  implementation  is  in  seven  files.  These  are  given  in  the  usual 

Pascal  ordering,  where  entities  are  defined  before  being  used  (or  referenced). 

The  main  program  is  contained  on  the  last  file. 

Your  attention  is  directed  to  an  earlier  section  (4.7),  where  it  is 
stated  that  this  implementation  will  not  provide  correct  results  for  touching 
flames.  Further,  the  default  values  for  the  model  parameters,  contained 

herein,  are  outdated. 
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alpha  :=  acos ( (SQR(rl)  -  SQR(d)  -  SQR(rn) )/(2.0*d*rn) )  ; 
xi 2  :=  xn  +  rn  *  COS ( theta+alpha) ; 

yi2  :=  yn  +  rn  *  SIN  (theta+alpha) ; 

xil  :=  xn  +  rn  *  COS (theta-alpha) ; 
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END; 
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•  begin 

promptat (quesline, 'ydistO  =  ?')#  READLN (ydistO) ; 
IF  ydistO  >  ydistf  THEN  ydistf  :=  ydistO; 
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promptat (quesline, ' azimuthO  =  ?  )»  READLN (thetaO) ; 
IF  thetaO  >  thetaf  THEN  thetaf  :=  thetaO; 
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DISTRIBUTION  LIST 


DEPARTMENT  OF  OEFENSE 

DEPARTMENT  OF  THE  ARMY  (Continued) 

Defense  Intelligence  Agency 

US  Army  White  Sands  Missile  Range 

ATTN 

RTS-2A,  Tech  Lib 

ATTN:  STEWS-TE-N 

ATTN 

RTS-2B 

USA  Missile  Comnand 

Defense  Nuclear  Agency 

ATTN:  DRSMI 

ATTN 

NAWE 

ATTN 

STSP 

DEPARTMENT  OF  THE  NAVY 

4  cys  ATTN 

STTI-CA 

David  Taylor  Naval  Ship  R  4  0  Center 

Defense  Technical  Information  Center 

ATTN:  Tech  Info  Ctr,  Code  522.1 

12  cys  ATTN 

DD 

Naval  Research  Laboratory 

Field  Comnand,  Defense  Nuclear  Agency 

ATTN:  Code  2627,  Tech  Lib 

ATTN 

FCPR 

ATTN 

FCT 

Naval  Surface  Weapons  Center 

ATTN 

FCTEI 

ATTN:  Code  X211,  Tech  Lib 

ATTN 

FCTT 

ATTN 

FCTT,  W.  Summa 

Naval  Surface  Weapons  Center 

ATTN 

FCTXE 

ATTN:  Tech  Library  4  Info  Svcs  Br 

Under  Secy  of  Defense  for  Rsch  4  Engrg 

Naval  Weapons  Center 

ATTN 

Engr  Tech,  J.  Persh 

ATTN:  Code  266,  C.  Austin 

ATTN:  Code  3263,  J.  Bowen 

DEPARTMENT  OF  THE  ARMY 

ATTN:  Code  343,  FKA6A2 ,  Tech  Svcs 

BMD  Advanced  Technology  Center 

Naval  Weapons  Evaluation  Facility 

ATTN 

ATC-0,  F.  Hoke 

ATTN:  Code  10,  Tech  Lib 

ATTN 

BMDATC-D,  M.  Capps 

ATTN 

1CRDABH-X 

Space  4  Naval  Warfare  Systems  Comnand 

ATTN:  PME  117-21 

BMD  Program  Office 

ATTN 

DACS-BMT 

Strategic  Systems  Programs,  PM-1 

ATTN 

DACS-BMZ 

ATTN:  NSP-43 ,  Tech  Lib 

BMD  Systems 

Comnand 

DEPARTMENT  OF  THE  AIR  FORCE 

ATTN 

BMDSC-HW 

Aeronautical  Systems  Division 

Harry  Oiamond  Laboratories 

ATTN:  ASD/ENSSA 

ATTN 

OELHD-DTSO,  00103 

ATTN 

DELHO-NW-P 

Air  Force  Institute  of  Technology 

ATTN 

OELHO-NW-RA,  L.  Belliveau,  22100 

ATTN:  Library 

ATTN 

DELHD-TA-L,  81100,  Tech  Lib 

Air  Force  Weapons  Laboratory 

US  Army  Ballistic  Research  Lab 

ATTN:  NTE 

ATTN 

DRDAR-BLA-S,  Tech  Lib 

ATTN 

OROAR-BLT 

Foreign  Technology  Division 

ATTN:  NIIS  Library 

US  Army  Engr  Waterways  Exper  Station 

ATTN 

Library 

Strategic  Air  Comnand 

ATTN 

WESSE 

ATTN:  DOCSD 

ATTN:  DOWE 

US  Army  Foreign  Science  4  Tech  Ctr 

ATTN:  INAO 

ATTN 

DRXST-SO 

ATTN:  NRI/STINFO 

US  Army  Material  Command 

ATTN:  DRXAM-TL,  Tech  Lib 

DEPARTMENT  OF  ENERGY 

Department  of  Energy 

US  Army  Mobility  Equip  R&D  Comnand 

Albuquerque  Operations  Office 

ATTN 

DRCME-WC,  Tech  Lib 

ATTN:  Technical  Library 

US  Army  Nuclear  4  Chemical  Agency 
ATTN:  Library 

US  Army  Tank  Automotive  R4D  Comnand 
ATTN:  DRDTA-UL,  Teen  Lib 


University  of  California 
Lawrence  Livermore  National  Lab 

ATTN:  Technical  Info  Dept  Library 

Oak  Ridge  National  Laboratory 

ATTN:  Central  Rsch  Library 


